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THE APPLICATION OF SPREADING 
MEASUREMENTS ON OIL TESTS. 


By N. E. M. Hacernorn and F. H. Stiettsss.* 


Following a description of the method of testing the spreading value of oil, 
a theoretical conclusion based on experiments made proves that the measure- 
ments as such are justified and, finally, the possibilities of putting this method 
into practice are discussed. 


INTRODUCTION. 


Ir a pure liquid non-volatile hydrocarbon is placed on watef, this will 
remain there in the form of adrop. A pure mineral oil acts the same way. 
Should the oil, however, contain compounds which are partly hydrophile 
and this oil is spread out on water, then these molecules tend to accumulate 
in the interface while they adjust themselves in such a way that the hydro- 
philic part of the molecule comes into contact with the water and the 
hydrophobe part turns away from the water. In this manner an oil film 
comes into existence. This film is not monomolecular, one gets a rather 
thin layer of hydrocarbon almost devoid of hydrophilic substances in contact 
with the water via the interface described above. This is called by Lang- 
muir a duplex film. The expanse of the surface of this oil film is a measure 
for the quantity of these compounds which are present. This is expressed 
in the spreading constant, which shows the surface occupied per unit of 
weight at the moment when the film just has started to exert pressure. 


MEASURING APPARATUS AND METHOD OF MEASUREMENT. 


The measurements are executed by the apparatus constructed by 
Gorter and Seeder. ! 

A barrier floats in a rectangular paraffinated glass container filled with 
0-001n-hydrochloric acid, to which container it is attached by means of 
very thin strips of platinum, so that it can move freely within certain 
limits. Forces affecting the barrier are transferred to a torsion balance 
and compensated by the torsion of a spring, care being taken to see that the 
barrier always returns to zero position. The amount of torsion is regis- 
tered on a scale, thus measuring the pressure on the barrier. Displacement 
of the barrier is indicated by means of an optical system. 

The oil to be measured is dissolved in petroleum ether and blown on the 
water surface between a paraffinated glass slide and the barrier by means 
of a micro-pipette. The glass slide is driven in the direction of the barrier, 
and as soon as the film starts to exert a pressure on the barrier, the torsion 
balance is adjusted as described above and the pressure, as well as the 
distance from the glass slide to the barrier, is noted. After this, the glass 
slide is moved closer to the barrier and the measuring process repeated. 
In this way a pressure curve is obtained. In most instances—between 
wide limits—this is a straight line. By continuing this line one extrapo- 
ates to pressure zero. The surface thus found is used for the determination 
of the specific spreading value. 

* Laboratory of the Netherlands Cable Works. 
UU 
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The petroleum ether used should be purified to such an extent that it 
does not affect the readings. This is obtained by repeated distillation, 
At the same time the water surface should likewise be so purified that no 
interfering molecules are present. If the spreading capacity of the oil is 
small, 0-05 ml petroleum ether is brought on the surface, which latter 
contains +15 mg of the oil. For oil which has a bigger spreading action 
a lesser quantity is sufficient. The determination takes three minutes, 
It is very fortunate that only small samples are needed. 


THEORETICAL CONSIDERATIONS. 


When making tests the following can be observed. After blowing the 
solution over the surface and moving the slide until the film exerts pressure, 
the major part of the surface is covered by the duplex film, there is, however, 
a water surface left uncovered. If the film is neither too thin nor too thick, 
an interference colour can be noted, indicating that the film is equally thick 
all over. After waiting a little, the picture changes. Holes appear in the 
film and the latter continues to contract, with local thickenings in it, often 
indicated by different interference colours. Provided one waits long 
enough, the film contracts to very small lenses. 

Considering the problem of the duplex film thermo-dynamically, one 
gets various conditions of equilibrium. One is mechanical and two others 
physico-chemical. The mechanical equilibrium of the surface tensions 
and hydrostatic pressures adjust themselves immediately, unless strongly 
viscous liquids are concerned. In the case of solids a mechanical equilibrium 
is also quickly established, which equilibrium may include anisotropic 
tensions, which makes it difficult to make any reliable calculations. The 
first physico-chemical equilibrium condition stipulates a certain relationship 
between volume concentration of the substances with hydrophilic groups 
in the oil layer and the surface concentration of these substances in the 
interface oil-water. According to Langmuir this equilibrium is established 
quickly enough, which is confirmed by our own experiments. This can 
easily be understood when considering that diffusion only takes place 
across the thin oil layer. The second physico-chemical equilibrium 
condition establishes a relationship between the surface concentrations of 
the active substances in the interface air-water and oil—-water. This 
equilibrium adjusts itself much more slowly than the other since here it is 
a question of diffusions along the surface. Langmuir? covers this in his 
experiments with stearic acid in petroleum. Our own experiments confirm 
this also, as described above. 


Mechanical Equilibrium. 


In the middle of the oil lense, where the interfaces show no curvature and 
are horizontal, the following horizontal forces work in the oil : 


1. The surface tension of the surface oil-air y,, (a tensile force in 
dyne/cm) ; 

2. That of the interface oil—-water y,,. (likewise a tensile force in 
dyne/cm). 
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3. The integrated lateral hydrostatic pressure, being equal at the 
top to the atmospheric pressure and increasing with gp, x dyne/cm? 
as one descends x cm (g is the acceleration of the force of gravity in 
cm/sec?, gg is the density of the oil). 


If there is equilibrium, and if the oil lense is large enough, these tensions 
must be equal to those in the adjoining oil-free water (the surface of which, 
however, may possibly be occupied by active substances) at the point 
where the water-air interface is horizontal. At the same time they must 
also equal those in the pure water at the other side of the barrier plus the 
pressure exerted per cm by the barrier itself. This results in : 


PolPw — Po) 


° P= Yua = (Yw)p —F.. (1) 


Yow + You ~~ hg 
t = thickness of oil lense in em; 
Yu)p = surface tension of the pure water on the other side of the barrier 
in dyne/cm ; 
ua = Surface tension of water—air interface adjoining the oil in dyne/cm ; 
F = the surface pressure exerted by the barrier in dyne/cm. 


Yow 80 Yyq Will depend on the concentrations of the active substances in 
the surface concerned ; +, does not alter, according to Langmuir,® when the 
volume concentrations of the active substances in the oil vary. The 
thicknesses of the oil layer with which we are dealing run from 0-1 » to 


several ’s, so the term : 
jg P0lPe — Po) yo 
Pw 


which is the difference between two integrated hydrostatic pressures 


Igpof® and does (£1) 
is approximately equal to : 
= J. 1000. 0-1 (10-5 till 10-3)? = 4(10-8 till 10-4) dyne/em 


so that this term can very often be disregarded. 

The above-mentioned thicknesses can be deducted either from the 
interference colours or by means of a simple calculation—i.e., dividing the 
total oil volume by the surface occupied. 

If the term is disregarded, it can be assumed that irrespective of the 
thickness of the oil-layer, the following equation applies : 


2 
’ 


Second Physico-chemical Equilibrium. 


At a certain barrier pressure F’, y_ the surface tension of the contami- 
nated surface water-air is fixed (see equation 1). Thus also the concen- 
tration of active substances on this surface. If the second physico-chemical 
equilibrium holds, the concentration of active substances in the interface 
oil-water should also be fixed, and therefore also the interface tension and 
at the same time the thickness of the oil-layer. 
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This conflicts with observations made. For instance, if oxidized oil 
is diluted with the original non-oxidized oil and this mixture is spread, 
entirely different thicknesses of the oil layer are obtained at the same barrier 
surface pressure than if undiluted oxidized oil is spread. 

It was also found, as has already been stated, that, provided one waits 
long enough, holes appear in the lense, the oil itself also changing in shape, 

The interesting thing is that in all probability it is the peculiarity of the 
films of not attaining the second equilibrium during the measurements that 
renders it possible to obtain simple and comprehensible results. It is 
essential, however, that measurements are made quickly. 


First Physico-chemical Equilibrium. 


Using the described method of measurement, there is always a rather 
small water-air surface present between the barrier and the glass slide. 
The major section of the surface is, however, covered by the spread oil 
film. If, therefore, the water-air surfaces are disregarded, and it is as- 
sumed that there is a physico-chemical equilibrium between the surface 
concentration of the active substances in the interface oil-water and the 
volume concentration in the oil, one can easily explain theoretically the 
relation revealed by experiments between : 


1. The surface occupied at certain barrier pressures (including that 
obtained by extrapolating to zero pressure) ; 

2. Total quantity of oil spread ; 

3. Quantity of active substances in the oil. 


This naturally gives a great deal of support to the acceptability of the 
suppositions. 

At a certain barrier pressure F there should be a certain y,,, (equation 
2). This means that if there is only one active compound present, the 
concentration from this compound in the interface oil-water o,, (molecules 
per cm?) is fixed. 

Similarly the volume concentration u, (molecules per ml) is also fixed. 
But then the following formula holds independent of the original concen- 
tration of active substances in the oil. 


A = dS + Uy V _. + + >» » 
in which we have constants o,,, and %, when the barrier pressure is fixed. 


A = the total number of active molecules in the oil; 
S = the surface in cm? of the oil; 
V = the volume in ml of the oil. 


From formula 3 it is easy to explain the seemingly paradoxical experimental 
results that, by adding pure oil (increasing V) to a certain quantity of impure 
oil (A therefore remains constant), a smaller surface is obtained with the 
same barrier pressure. This means, therefore, that the thickness of the oil 
layer increases out of proportion to the volume. The new pure oil added 
must also be brought up to the concentration u, of active molecules. These 
active molecules can only be obtained from the interface, which must, 
however, retain the same concentration, so that the surface must of neces- 
sity become smaller. 
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Before we discuss the results of our experiments any further, we amend 
the formula (3) somewhat and go on to specific quantities : 


S_ 1A 
a a 


> is the specific spreading ;: Surface in cm? per total volume of oil in ml. 


We will call this s (usually the specific spreading is given in cm?/mg oil, 
but for a theoretical research it is better to use the cipher cm? per ml). 


= is the total concentration of active molecules per unit of volume (thus 


including those present in the dividing surface). This number is unknown, 

and therefore we substitute a dilution z. 

We started with a certain oxidized oil, which we assume had a total 
concentration a, of active molecules per ml. 

Spreading tests were made with this and other oils, which were obtained 
by mixing the latter with varying quantities of the original pure oil con- 
taining no active substances. 
oxidized quantity of oil 

total quantity of oil 
Formula (4) thus becomes : 


The relation we will call dilution z. 





Fig. 1 shows s as a function of x. There are three curves. The first 
stands for the extrapolated barrier pressure of 0 dyne/em; the second for 
a barrier pressure of 1 dyne/cm and the third for 1-5 dyne/em. It can be 
seen that the points fall quite nicely along straight lines. If these lines 
are continued, they intersect the s axis. —_ section between the 


inter-section point and the origin should equal —. The reciprocal is a 


linear dimension in cm. This magnitude equals “the thickness of an oil 
layer containing as many molecules as the corresponding interface oil- 
water. Corresponding means that the concentrations at the interface and 
in the oil are in physico-chemical equilibrium with each other at the given 
barrier pressure. In the case in question this thickness appears to be 
~10-4em = In. 

When x = 1 (i.e., the unmixed oxidized oil), the actual thickness of the 
oil layer at a pressure of zero is ~0-06.; the major part of the active mole- 
cules is therefore present in the interface. 

It is interesting to note that the lines for 0, 1 and 1-5 dyne/cm all intersect 
the s axis at the same point. This means that at the equilibrium there is a 
constant ratio between the volume and surface concentrations for these 
pressures, which gives strong support to the hypothesis that the active 
molecules in the interface are present in a two-dimensional gas state in 
equilibrium with non-associated molecules in solution in the oil. This 
hypothesis complies with Gibb’s theorem, and demands a conforming with 
Yoa + Yow — Yu, Of the pure oil in the same manner as Langmuir * 
has done with stearic acid and oil. It should be remarked that Langmuir 
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gets quantitive results different from that which we find with oxidized oils, 
This leads us to believe that the active molecules, which cause the previously 
mentioned relations, are not acids. This is also apparent from the small 
influence of the pH of the water shown in the result of our tests. Should 
such acids still be present, however, it can be understood that they would 
be in a condensed two-dimensional phase, and would exert a small, constant 
saturation pressure. 

Another consequence of the constant ratio between the surface- and 
volume-concentration can be deduced from formula (3). 





s incm2?/cm3 
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the surface axis to smaller surfaces when the oil is more diluted. 
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This rather remarkable phenomenon was the first indication we received 
of the existence of the simple relationships described above. 

If one examines what happens when two kinds of molecules with different 
surfaces are present, while the two-dimensional gas law applies, and at 
the same time these molecules are dissolved in non-associated condition, 
one obtains that the resultant theoretical relation between s and x will not 
be rendered by a straight line, but by a hyperbola. Over a great stretch 
this hyperbola cannot be distinguished from its asymptote, so that the 
straight line connection can still be used. (The straight line applies when 
the oil-layer is so thin that the major part of the active molecules are found 
inthe surface.) The theoretical curve (the other asymptote of the hyperbola 
runs parallel to the x axis at a certain distance below it) likewise dips below 
the « axis, where x has a positive value. In fact, of course, this is not 
feasible. The digressions of the actual curves as compared with the 
theoretical ones are caused by the fact that in the case of excessive dilutions : 


(1) That part of the space between the barrier and the strip occupied 
by the free water surface can no longer be disregarded ; 

(2) The oil layer becomes so thick that the simplified equation (2) 
is no longer applicable. 


In any case one never gets too high values when the specific spreading 
(surface at barrier pressure zero divided by the total quantity of oil) is 
taken as a measure for the content of active substances. This content is 
always bigger in reality and the error relatively smaller as the content 
increases. 


APPLICATION. 
Cable Research. 


In such research the spreading values often give good indication as to 
what is happening in the insulation. The insulation of a high-tension cable 
consists of paper with impregnated oil. The oil used for this purpose must 
be of very special quality, and it is obvious that very thorough methods of 
investigation must be used in order to trace any changes in the oil. 

The curves a of Fig. 2 show clearly that it is possible to point out very 
small differences. We are dealing here with centrifuged oil from the 
paper from different layers in the insulation of four pieces of cable. Of 
these, III and IV have been kept at an increased temperature for a certain 
length of time, while I and II have remained at room temperature for the 
same period. Furthermore, II and IV have been subjected to electrical 
tension tests. The spreading values are shown as a function of the number 
of the layer. 

In addition to the spreading measurements, the power factor of the 
oil-impregnated paper is measured with a condenser such as described by 
Wyatt.4 In curves b of Fig. 2 these measurements are shown. These 
paper tapes, after being measured at 75° C, were centrifuged too and the 
spreading values of the oil measured. See curves c of Fig. 2. 

The remarkable thing is that the power factor measurements do not 
show the differences between the four pieces of cable in the middle of the 
insulation, found with the spreading measurements. 
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Fie. 2. 


MEASUREMENTS ON FOUR DIFFERENT CABLES PLOTTED AGAINST THE NUMBER 
OF THE PAPER LAYERS COUNTED FROM THE LEAD SHEATH. 


Cable kept at room temperature. No electrical tension applied. 
-<------ Cable kept at room temperature. Electrical tension applied 90-100 kV. 
—+— +. Cable kept at 55-75°C. No electrical tension applied. 
tee eeeeee Cable kept at 55-75° C. Electrical tension applied 90-100 kV. 
a: Specific spreading values of the oil centrifuged from the paper tapes. 
b: Power factor measurements of the impregnated paper tapes at 75° C and a 
frequency of 76 c.p.s. 
C : Specific spreading values of the oil centrifuged from the paper tapes after power 
factor measurements. 
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In addition, it is also noticeable that the oil in the condenser oxidizes 
slightly at a temperature of 75° C. 


Oil Oxidation. 


When making oxidation tests we use the specific spreading as a measure 
for the grade of oxidation. The top graph of Fig. 3 shows such a process 
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Fie. 3. 


SPREADING VALUES AND OXYGEN ABSORBED AS A FUNCTION OF THE TIME IN 
HOURS, 


with oxygen absorption plotted against the time. The way in which these 
oxidation tests are made has, of course, no bearing on the subject. Our 
test can determine the exact moment at which the oil commences to absorb 
oxygen, and in this particular instance the process began 15 minutes after 
the start of the test. It is interesting to know if anything can be noticed 
in the oil in this preceding induction time. Spreading measurements give 
an affirmative answer. When the process began, the values were measured 
every ten minutes; they are reproduced in the lower graph. It can be 
observed that the oil changes already in the induction time. Now compare 
these spreading measurements with the acid value. After an oxidation 
time of two hours the oil has a specific spreading value of 20 cm*/mg, the 
original oil having a zero value. This is certainly no small value for the 
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oxidized oil. The acid value is found to be 0-07 and 0-05 mg KOH /g oil 
for the oxidized and non-oxidized oil respectively. This result is self. 
explanatory. 


Purifying of Oil. 

It is possible to remove oxidation products from oil by treating it with 
charcoal or clay. The effect of such a purifying process can easily be 
followed by the aid of the spreading.5 

Oil with a spreading value of 52 cm*/mg was treated several times 
consecutively with charcoal and clay (10 per cent by weight) at a tempera- 
ture of 100° C. In the case of charcoal, after one treatment the value 
dropped from 52 to 33, after the second treatment to 7, and after the third 
to 3cm?/mg. With clay the results were resp. 4, 2, and 1 em?/mg. ‘The 
Waterman ° analysis shows that the composition of the oil changes less 
when treated with clay than when treated with charcoal. 


Examination of Hydrogenated Oils. 


Various samples of fractionated distillates were measured as well as the 
products after hydrogenation. The original products had a positive spread- 
ing value, which increased, according to the length of time the fractions 
had been kept. When hydrogenated, the active molecules disappear. 
In tests made, complete disappearance took place, however, only after the 
oil was entirely hydrogenated. In case of partial hydrogenation some 
active molecules are still present. 

Here follow two examples : 

First a fraction with its hydrogenated products, after having stood for 
a period of six years. The original oil had a value of 755 cm?/mg; the 
partially hydrogenated a value of 360 cm?/mg and that entirely hydro- 
genated a value of 3 cm?/mg. 

In comparison to this a measurement was made of a fraction and its 
entirely hydrogenated product immediately after distillation. The values 
were respectively 14 and 0 cm?/mg. This clearly illustrates the influence 
exerted by time. 


- 


CONCLUSION. 


The spreading measurements described above offer a particularly sen- 
sitive method for investigating the properties of mineral oils which are of 
importance for practical use from the standpoint of their oxidation stability, 
purity, etc. This method often enables one to show changes where other 
methods fail. The theoretical considerations teaches us to understand 
some of the formulae reached, and arouses the expectation that the method 
could be still further developed along the lines of what might be called 
differential diagnoses. 

Finally we wish to pay a tribute to the memory of the late Mr. H. W. G. 
Van Gils, E.E., who instituted this method of research in connection with 
insulating oil at the laboratory. He was killed as a civilian by the 
Germans on May 10, 1940. 

Due to the isolation caused by the war, we only learned of the papers 
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by W. A. Zisman and by W. D. Hartins (J. Chem. Phys., 1941) after the 
above paper had been written. We think that the interesting general 
investigations in these papers are supplemented from the practical 
viewpoint by our research. 
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FRACTIONAL DISTILLATION OF TERNARY 
MIXTURES. PART II. 


By A. J. V. UnpERwoop, D.Sc., M.I.Chem.E., A.M.I.Mech.E., F.R.I.C., 
F.Inst.F. (Fellow). 


SuMMARY. 


The basic equations derived in Part I are employed to determine the con- 
ditions for a state of minimum reflux. Equations are presented by which 
the minimum reflux ratio for ternary mixtures can be readily calculated when 
the fractionation between the key components is a sharp one. 


In Part I of this paper} equations were presented for computing the 
composition on any plate in a fractionating column for ternary mixtures, 
Since it appeared an alternative method of making these computations has 
been described by Harbert.? 

An alternative, and rather more elegant, derivation of the basic equations 
presented in Part I is the following. 

The compositions of the liquids on adjacent plates are given by the 
relations 


ee re || 
_—* yX, + BY, + % 


o a 
mye tem + By, +2 


(3) 


As previously, for a rectifying column, m = rev where FR is the reflux 
a a -_ 2 _ 
Rol C=Ra1) d Rai1°” and 6 are the relative 
volatilities of components x and y to component zandy >8> 1. Fora 

S+1 


stripping column, m = = where § is the “ reboil ratio 


number of moles of vapour produced in the reboiler per mole of bottom 


ratio; 6 = 


, 


’—4.e., the 


product withdrawn; 6 = ; = 3 — =. For a stripping 


column on which a rectifying column is superimposed, S = 


For both rectifying and stripping columns 
b+c+d=1—~Mm. 


Constant molal reflux and constant relative volatilities are assumed. 
— ee 
y—-¢ B-¢ 


Now multiply equations (1), (2), and (3) by and - 





m 
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respectively (where ¢ is a quantity as yet undetermined), and add. 
Then 

Bo 4% 
ol tera tes 





Now choose ¢ so that 








y%o BYo 
y—9" F-8" 
Equation (4) is of the third degree, and gives three values of ¢, denoted by 
$1, $a, $3, it being understood that ¢, < ¢_ < $3. 

Substituting these values in equation (5) and applying the method used 
in Part I, there are obtained, for the composition on the nth plate, the three 
equations 





= ae myx, + By, + 2) 


Y%o_ 4 _BYo_ ¥%n a 
rey ha cy ba _(#" y-& 8-6 1-4 
_ Xo _ BY Yrn BYn zn 
y—$ | B—ds y-h 8-4 1-4 
. Y7o __ BY 20 YXn BYn 
1-4 g— 4, 71-4 _ (hy yw Bw” 
Yo _ + BY + 20 3 Y=n + BYn 
y—¢ B-¢ 1-4 y—¢ B-—ds 
Yo _ BYo 29 Y2n BYn Zn 

rr Ae eT & _(#) Yoh B-B IKK @ 
Yo, _BY 2 , 4 Pte 
y—-¢, B-—4¢ — $y y—-¢, B-¢ 1-4, 


In the above derivation the only property assumed for the parameter ¢ 
is that it should satisfy equation (4). The other properties of this para- 
meter, stated in Part I, can be readily deduced. Assume the composition 
of the liquids on any two adjacent plates in the column to be the same, as 
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is the case under minimum reflux conditions. If this composition js 


denoted by (h, k, l), then, from equation (5) 
@$=mMyh+BkEK+1)... . =. =. (7) 
Equations (1), (2) and (3) then give 
mh+b mk+c_ mi+d _ 1 _m 
yh BR i nh + BE +L PF 


a ae ee 
= nly — 8) mip — 9) ~mi—g 
For the three values of ¢—i.e., 6), $2, ¢;—Wwhich satisfy equation (4) there 
are three values of (h, k, 1) which will be denoted by (h,, k, 1,), (hg, ke, Is) 
and (hg, kg, /,), respectively. 
For a rectifying column 6, c, d are positive, and it has been shown in 
Part I that, in this case, 


0<¢,<1; 1<4,<8; B<¢,<y. 


For ¢,, equation (8) will give positive values for h,, k,, and l,. For ¢,, 
h, and k, will be positive and /, will be negative. For ¢5, A, will be positive 
and k, and /, will be negative. 

For a stripping column, it will be convenient to write 5, ¢, d instead of b, 
c, d and ¢ instead of ¢ and h, &, ] instead of h, k,l. In this case §, ¢, d are 
negative and it can be readily shown that 


Ll<$ <8; B<d<y; ¥ < 4. 


For },, equation (8) will give }, and &, negative and ], positive. For 4,, 
he Will be negative and {, and /, positive. For 4s, hs, ks, and ], will all be 
positive. 

For both a rectifying column and a stripping column there are thus three 
cases where compositions on adjacent plates are the same, and each of these 
cases will correspond to minimum reflux conditions. Equation (4), 
applied to a rectifying column, can be written in the form 








yD BYD +, | i 


T=9. §—@ a 
For any given value of R, equation (9) gives the corresponding value of ¢, 
and those of h, k, 1 are then obtained from equation (8). These values of 
h, k, | represent the limiting compositions which can be attained with the 
given reflux ratio when the number of plates is infinite. 

The various types of cases which arise can be illustrated by the following 
three examples which all deal with a rectifying column. In all three 
examples y = 4, 8 = 2, R = 3. 

Example 1. All three components are present in substantial amounts in 
the top product. xz» = 0-4; yp = 0-4; zp =0-2. Equation (9) gives 
¢, = 0-9267; ¢, = 1-774; ¢, = 3-649. Equation (8) gives 

1 = 0-040; kj =0-115; 1, = 0-843 
2 = 0-106; kz, = 1-046; 1, = —0-153 
1-386; k, = —0-295; 1, = —0-092. 
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Example 2. The heaviest component is present in only a small amount 
in the top product. 2p = 0-599; yp = 0-400; zp = 0-001. 
Then 
900583 ; $2 = = ]- 7295 $y = 3-473 
k, = 0- 850; 1, = —0-000791 
k, = —0-3145; 1, = —0-000468 


Example 3. The two heavier components are present in only small 
amounts in the top product. 2p = 0-999; yp = 0-001; zp = 0-000001 = 
1 x 10°. An extremely small amount of the heaviest component is here 
included. In any fractionation of a ternary mixture all three components 
will be present in the top product. Where the second heaviest component 
(y) is present in small amount in the product, the amount of the heaviest 
component in the product is quite negligible for all practical purposes and 
is not normally taken into account. The effect of including it in the 
calculation will be discussed later. 

Then 


5 by = 1999; $5 = 3-0015 
k, = 0-00033; 1, = 0-88867 
k, = 0-667 ; l, = —0-67 x 10°6 
ks — —0-001; L, = —05 x 104 
The results for Examples 2 and 3 have been calculated to a degree of 
accuracy which is obviously greater than would be required in any practical 
problem. The many decimal places shown in some of the results do not, 
however, mean that the calculations require an accuracy in computation 
greater than that of a slide-rule. When any of the components are present 
in the product in small amount the calculations are considerably simplified. 
For instance, in Example 2, zp is small. The root ¢, of equation (9) is 
therefore very nearly equal to 1. Equations (8) can be written in the form 


oi. to, ee: ¢ . 2p 
*=Ry—@) “—Re—-9 '~RI-a °° 


As $, is very nearly equal to 1, the values of h, and k, become Ry =D and 


ip Since h, + k, +1,= 1, the value of J, can be found by 
difference. 

Alternatively, putting ¢, = 1 in the first two terms of equation (9), we 
have 


This gives the value of the small quantity (1 — ¢,) and also the value of 
i 5 which cart then be used to find J, from equation (10). For this 
PE 
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example the other two roots ¢, and ¢, are found by neglecting the term 
rs in equation (9) and solving the resulting quadratic equation. 

In Example 3 both yp and zp are small. In this case ¢, is very nearly 
equal to 1 and ¢, is very nearly equal to 8. For finding 4, and h,, k,, L 
the same procedure is used as for Example 2. For ¢,a second approximation 
is made as follows. Since ¢, is very nearly equal to 8, the term r “D r 
can be neglected, and equation (9) can be written (with 6 for ¢, in the first 
term) as 








YXp BYp 
=R 8 
eab ta + 


This gives the small quantity (8 — ¢,) and the value of P_ required 


2 
for finding k, from equation (10). It has been shown in Part I of this 
paper that ¢,¢.4¢, = m$y, and when ¢, = | and ¢,= 8 approximately, 

Y 
R+1 
can be found by a similar procedure to that already described for ¢, by 
writing 


then ¢, = my = 





approximately. A second approximation for ¢, 


Most discussions of calculations for minimum reflux conditions have 
dealt with the case of a rectifying column superimposed on a stripping 
column. This is obviously the most important type of case in practice, but 
certain important principles can be illustrated more clearly by considering 
a rectifying column or a stripping column alone. For the following 
discussion a rectifying column alone is considered, and is assumed to be 
mounted on a still kettle containing liquid of composition (xp, yp, z,). 
We now have to find the limiting kettle composition for a given reflux ratio 
when the number of plates is infinite. 

A method which has been used by Colburn,’ by Gilliland * and by the 
author > for the case where the top product contains one component in 
small amount is the following. If the compositions on two adjacent plates 
are the same, equations (2) and (3) give 





mr +e 8 YP 

mip +d" Zp 
‘ ‘ vr a ee: a - = o-— = 
and if d is negligible, then — = 1) ~* R@ —1) Similarly 
tp = By This method obviously gives the values (h,, k,, 1;) 


obtained by putting ¢, = 1 in equations (10). 

If yp is small as well as zp, the usual procedure has been to neglect zp 
altogether—a very natural thing to do, as zp must be very small indeed, 
and, if a small value yp is specified in the product, the value of zp cannot 














: 
t 
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also -be specified ab initio. Neglecting ,zp and also neglecting c in the 
equation 

mip +b _ yap 

myp +e BYP 





we obtain 


b Xp 1 


aE 1) R(Z—-1) RE ) 


(It is obvious that this represents the limiting composition under minimum 
reflux conditions for the binary mixture of components x and y.) It is 
then assumed that, with this value of xp, a “ pinched-in region ” occurs in 
which zp is negligible and yy = 1 — xp. The values so obtained are those 
of (hg, kg, 1) obtained by putting ¢. = $ in equations (10) and neglecting the 
small value of /,. The calculation can be carried farther by adding in a 
small amount of component z and carrying out a plate-to-plate calculation 
down the column. 

It would, however, be equally legitimate to base the calculation on zp 
instead of yp, and there would then be obtained the limiting compositions 
(hy, &,, 4) corresponding to ¢,. As will be seen from the results of 
Example 3, where zp = 1 x 10°, the final effect of including it in the 
calculation is very marked. The values of h,, ky, l, calculated for this 
example can, of course, also be arrived at by making a plate-to-plate 
calculation down the column, including the component z in the calculation. 

The method just described does not give a correct conception of minimum 
reflux. It fails to provide answers to the following questions :— 

(1) In the apparently straightforward case where only zp, is small the 

<p Yb 
~ Ry —1) “= RG 1) 
which is independent of R. It therefore 





limiting compositions are given by xp = 


xp(8 — 1) 

ynly — 1) 
depends only on the a. of the product and the relative volatilities 
of components x and: y, and is the same for any reflux ratio. If the still 
kettle contains a charge in which the ratio of these components is other than 
r,(8 — 1) 
Ydly te 1)’ 

(2) If both yp and z, are small, very different results are obtained for the 
limiting compositions, according as zp is omitted from or included in the 
calculation. Which of these results is correct ? 

It is believed that the following analysis clarifies the position. 

It has been shown in Part I that if we put 29, yo, z) equal to rp, yp, Zp, 
respectively, the left-hand side of equations (6a), (6b) and (6c) becomes 
unity. From equations (6a) and (6c) the composition on the nth plate 
below the top of the column is given by 


The ratio of xp to yp is 


how does this method of calculation apply ? 








YXn + BY, an 
y—¢. | ey ed 1 —¢. __ ¢,\" 

Yen, _BYn 4 2 nt: (31) » «+ « (Ila) 
Y¥-@ B6- $y Fe a 
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(12a) 





(12b) 


It is necessary that the values of 2,, y,, 2, which satisfy equations (12a) 


and (126) should not make the denominator sere + 3 rie + ra ; 
equal to zero. The values of x,, y,, z, which meet these requirements on 
hy, ky, l,. 
The condition that h,, k,, 1, should satisfy equation (12a) is 
yhy Bk, l, 
+ - +——=|—~_=0... . 13 
y—h '8-& 1-4 i 


___ __ __ 4 
rom easton Oe = any = Ga) | mB — Hy wal =F) 








and 
equation (13) becomes 


cB d oe 
+ al — a) ~ 








0. (14) 


by + 
(y — daly — 4) (8B — $e)(8 — 4) 


Now, from equation (4) 


- an. See. by cB = 





fuk bh to eek eH 


by(di — 2), _ Bld — dy) d(¢,—%) _ 9 
(y — daly —$1) (8 —be)(B —4 3) (1 — $)(1 — 44) 


which is the same as equation (14), so that h,, k,, 1, satisfy equation (12a). 
In exactly the same way it can be shown that they satisfy equation (120). 
Also 


yhy + BK + hy sia bydy cBd, dy 
y—¢, B-—¢ I1-—¢, my—¢, m(B — ¢, m(1 — ¢,)* 


All three terms in this expression are positive for a rectifying column or 
negative for a stripping column, and therefore their sum cannot be zero. 

It can similarly be shown generally that the three equations of the type 
of (12a) are satisfied as follows :— 


a 








yt yt 





yx by er 
v4, 15-4 Ie TO 





\ 
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—~__ = 0 by hy, ky, 1, and hy, ky, ly 
— $y 


— 3 


If any two of the equations are satisfied simultaneously the solution is 
the values of h, k, 1, which correspond to the value of ¢, which does not 
appear in the two equations. 

Suppose now that a plate-to-plate calculation is made starting from the 
top and proceeding upwards. This can be done by putting n = —N (where 
N is positive) to give the composition on the nth plate above the top. From 
equation (6a), (6b), and (6c) the following equations can be derived in the 
same way as equations (lla) and (116). 








= 0 by hy, ky, ty, and hy, ky, le 


Yt-n BYy_y 
y— $) ui . ~ 


y*-N 





yn Syn ZN 
y—-& 8-& 1-4 _ 
yey Pw | an 


rN Ai qunnles 
yh | 8-& 1-4 
When N becomes infinite, 





ys -, Per . S77 9 |. . (Oe 
I-k be Toe — 


Yy*-n . by_y . — “iy 160 
y— Shah Soe = 


The solution of these equations is then (hg, ks, /,). 

Thus, if a plate-to-plate calculation is started at the top (including all 
components) and continued downwards, the limiting composition finally 
reached is (h,, k,, 1,). If the calculation is made upwards, the limiting 
composition finally reached is (hs, ks, Is). 

We now have to consider how the limiting composition (hg, ky, 1.) can be 
reached. Suppose now that the values of (x, y, z) on any plate are such 
that they satisfy the equation 





= by aa. . . (7 
y — 3 +ej-|GTI- ; (17) 


but do not satisfy the other two equations, so that these values of (x, y, z) 
do not represent limiting compositions. From equation (5) it is seen that, 
if equation (17) holds good for a plate it will also hold good for the next plate 
above or below it, and therefore for all plates above or below it. From 
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equation (6a) it will be seen that, as n increases, the left-hand side of the 
equation continuously decreases, and "a becomes zero when 1 is 


infinite. This requires = 4 PY $ ~ ag =0. As equation 
(17) is also satisfied, this means that the = ied composition (hg, k,, 1.) is 
reached as the calculation is made up the column. Note that hg, kg, 1, do not 
make the denominator of the left-hand side of equation (6a) zero. In 
exactly the same way it can be shown that the limiting composition 
(hy, ky, 1) is reached as the calculation is made down the column. 

Thus, if the composition of the charge in the still kettle is such that it 
satisfies equation (17), a plate-to-plate calculation up the column, starting 
from the kettle, will finally reach the limiting composition (hg, kg, 1,). A 
slight increase in the reflux ratio will result in an increase in the values of ¢ 
which satisfy equation (9), for, differentiating this equation, we have 


ie, 
orto =e tier *0 hts = - 








With an increase in the value of ¢, the expression — o $, + a are + i — J 
3 “oe — 
will no longer be zero, but will have a positive value. The plate-to-plate 
calculation up the column will then no longer pass into and remain stuck at 
the limiting composition (h,, kg, 1,), but will be capable of being continued 
until the composition corresponding to the top product is reached. Thus 
any composition of the kettle charge which satisfies the equation (17) 
will be a composition for which the reflux ratio used in calculating 
the values of ¢ from equation (9) is the minimum reflux ratio. The two 
limiting compositions (h,, k,, /,) and (hg, ky, ly) are merely particular values 
which satisfy the general equation (17). This conception of minimum 
reflux appears to answer satisfactorily the two questions which were stated 
earlier in this — 
Examples 1, 2, and 3 show that the value of bl is always negative. 
Although a slight increase in the reflux ratioto as ” Ae ra = J 
3 
positive will permit of a plate-to-plate aiabtien bole continued to the 
top of the column, it does not necessarily follow that this calculation will 
not pass through a negative value of z when plate compositions approaching 
hg, ke, ly are reached in the calculation. A plate-to-plate calculation which 
passes through a negative value of one of the components obviously does 
not represent a practical case. In a case like Example 1, where all three 
components are present in substantial amounts in the top product, /, has a 
substantial negative value, and when the reflux is increased slightly over 
the minimum, a plate-to-plate calculation could still pass through a 
negative value of z. This case is, however, an unusual one. In the normal 
cases, such as Examples 2 and 3, where one or two of the components are 
only present in small amount in the top product, /,, although negative, is 
very small. A quite minute increase in the reflux ratio above the minimum 
will ensure that the plate-to-plate calculation does not pass through a 
negative value. In other words, the theoretically correct minimum reflux 
ratio is that which will give a value of z (at the appropriate point in the 
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column) of —0-000791 in Example 2 or —0-67 x 10° in Example 3. The 
minimum reflux ratio for practical purposes is that which gives a value of 
zero instead. The difference is obviously outside the range of normal 
calculations. 

It is to be noted that, as in the case of a binary mixture, as long as the 
amount of a component in the product is small, the actual amount does not 
affect the minimum reflux ratio appreciably, as the limiting compositions 
do not occur in the neighbourhood of the product composition. If, however, 
a calculation is made for a finite number of plates, the actual amount of a 
component present in small quantity will appreciably affect the number of 
plates required for a given separation. 

It has been shown that, for a rectifying column, if the composition at any 


point in the column satisfies the equation i. BY aoe ae 


‘5-4 1-h 


there is a possible range of compositions between (hy, ky, ly) and (hg, kg, 1). 


It can be similarly shown that, if the equation - 2. + 
Y 


By Pie. 
. _ =0 
B—¢, ' 1—¢, 


is satisfied, there is a possible range of compositions between (h,, k,,1,) and 


yx By z 
_ and that if the equation —“— + —“— + -=3@ 
aaeolee ened WL 


is satisfied there is a possible range of poe ae between (he ky, 1,) and 
(hg, ky, ly). It appears, however, that these other two cases are not 
significant for practical problems of fractionation. 

The analysis has been given in detail for a rectifying column. It can be 
made in exactly the same way for a stripping column. For a stripping 
column, equation (4) becomes (writing Y instead of ¢) 


rate = 4 : 5: Lag va ae 
Here also }, < Yo < 4. 

For a stripping column, the limiting composition (h, ks, /;) corresponding 
to Y, is the one reached by calculating up the column from the reboiler, 
using all three components in the calculation. The limiting composition 
(hy, ky, 11) corresponding to , is the one reached by calculating downwards 
below the reboiler. The limiting composition (hg, kg, /,) is the intermediate 
one, and corresponds to that reached by calculation upwards from the re- 
boiler when two of the components are present in small amount in the 
reboiler and the lightest component is not taken into account in the 
calculation. 

For any given value of S in equation (18) minimum reflux conditions will 
obtain in the stripping column if the composition at any point satisfies the 
equation 





yx By z 
a — +: —=@O. 6 ow 19 
y—*4 B— 1—y, (19) 


and this equation represents a range of compositions between (hg, ks, /,) and 
(ha, ke, 12). Similar relations can be obtained for the other equations, but, 
as in the case of a rectifying column, they do not correspond to practical 
cases. 
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It should be noted that equation (17) for a rectifying column contains 
¢3, the value of which lies between 8 and y, while equation (19) for a stripping 
column contains },, the value of which lies between 1 and 8. 

For a stripping column, equations (8) become 


oon a a ees Lene ee. (20) 
(S+1)(y—¥)’ (S + 1)(8 — 9)’ (S+1)l—¥) 
For ordinary purposes it is not necessary to calculate compositions to 
several decimal places, as was done in Examples 2 and 3, and this results 
in considerable simplification. When x» = 0, equation (18) gives }, = y. 
Then, from equations (20), 


aE A, Ae. 
oS + Illy — 8) FS + Dy — DV) 
and, by difference, 
a: OO : 
hs S+)y—6) SO +Iy—) 
If, in addition, yw = 0, then equation (18) gives }, = 8 and h, = 0, 
eS Se. cee 
2 ($+ )(8—-1I)” S+1H8-—) 
since zy = l. 
: ee, Cae ee ee 
By difference, k= 1 +e —)’ 


If both zy and yyw are zero, equation (18) gives 


cw a ie S be ot. 
i. Sor }, = g ° 
and h, = 0, k, = 0, and 1, = 1. 

Similarly for a rectifying column, if zp = 0, equation (9) gives 4, = 1 
and equations (10) give 


a — ee 
= Ry —1y "= R@—N 
and, by difference, 


_ es. ee =e 

s Riy—1) R(®B —1) 
If, in addition, yp = 0, equation (9) gives ¢. = ® and, from equations (10), 
8 


l,=0O and, by difference, k, = 1 — ————. Also, 
. y differe ; Ry —8) 


6 
h, = =: 
2" Rly — 8) 
for this case, equation (9) gives 


yR 


= bs it =R a 1 and %=RiT 


i... 
Then &, = 1, &, = 0,1, = 0. 
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When a rectifying column is superimposed on a stripping column, true 
minimum reflux conditions obtain for the system when minimum reflux 
conditions obtain in both the rectifying column and the stripping column. 
With a mixture of three components, of which two are “ key ” components, 
one component will be small in the top product and two components will be 
small in the bottom product or vice versa. (The following analysis is 
limited to such sharp separations.) 

Consider the former case, in which the third component is lighter than the 
key components y and z. The top product contains component x and 
component y with a small amount of component z. For minimum reflux 
conditions in the rectifying column, by calculating downwards, the limiting 
composition h,, k,, 1, is reached. For the stripping column, the bottom 
product is component z with a negligible amount of x and a small amount 
of y. If on any plate below the feed-plate the equation 


tad By $s 0 
y—*" 5-a f-% 
is satisfied, minimum reflux conditions will obtain in the stripping column. 


If, at the feed plate, the values h,, k,, 1, for the rectifying column satisfy the 
equation 








2... Mh tet ent 


y—h B— vy 1—y 


for the stripping column, then minimum reflux conditions will obtain in 
both columns. The condition for minimum reflux is therefore 


Now it has been shown previously that h,, k,, 1, is a solution of two 
equations for the rectifying column—namely, 


yhy + Bk, l 
y—¢2 B — de — ds 


re ta Co he oo a 
There are thus three equations to be satisfied, and there are actually only 
two variables in these equations—namely, ? - t The condition that 
the three equations are satisfied simultaneously is that the determinant 


=0 








Ti 
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This determinant becomes zero if any two rows are the same, and this 
occurs if ), = $. or }, = ¢, or dg = $3. Now, ¢, cannot equal ¢, because 
¢, lies between 1 and § and ¢, lies between 8 andy. Also }, cannot equal 
¢3 because , lies between 1 and § and 4, lies between 8 and y. But 4, 
and ¢, both lie between | and 8, so that }, = ¢, is the only possible solution, 
This, therefore, is the condition for minimum reflux in the case where the 
third component is lighter than the key components. It is to be noted 
that ¢, is the root corresponding to the light key component in the rectifying 
column, the components being regarded in the order of increasing volatility 
and the values of ¢ in the order of increasing magnitude. In the same way, 
, is the root corresponding to the heavy key component in the stripping 
column. 

The other case to be considered is that in which the third component is 
heavier than the key components. In this case the key components are 
x and y. The top product is component x with a small amount of 
component y and a negligible amount of component z. The bottom 
product is components y and z, with a small amount of component x. For 
a stripping column the limiting composition, proceeding upwards, is reached 
at hg, ks, ls. For the rectifying column, minimum reflux conditions obtain 
if the equation 

~~ By ‘ 2 7 
y—¢ B-—¢s 1-45 
is satisfied, and this equation must be satisfied by h, ks, 1; for minimum 
reflux conditions to obtain in both sections of the column. Now, hg, ky) Is is 
given by the equations 








0 


By the same reasoning as before, the condition that all three equations are 
satisfied is }, = ¢,. Here again the value of ¢ corresponds to the light key 
component in the rectifying column and the value of 4 to the heavy key 
component in the stripping column. 

The general condition for both cases is that equation (9) for the rectifying 
column and equation (18) for the stripping column have a common root. 
Denoting this common root by 0, then 


yxp ByD 2 _ ip 
tet reaet atl...) am 





yew , BYw a 
Steet es 


Multiplying the first equation by P and the second equation by W and 
noting that 


Prp + Way = Fup; Pyp + Wyw = Fyp; Pzyp + Wew = Fry 
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where xp, Yp, Zp are the composition of the feed, then 


= +c + ar = (R + IP — SW 


Now, SW = RP+qF —W=(R+1)P—(l1—@)F 
since W = F — P. 
Equation (24) then becomes 


y<r BYP an 
es ak E51 all £21 Se 


1—6 


If the feed is liquid at boiling point, g = 1 and equation (25) reduces to a 
quadratic. If the feed is all vapour, q= 0. In this case equation (25) 
obviously has one solution @ = 0 since zp + yp + zp = 1 and equation (25) 
again reduces to a quadratic. 

If q = 1, equation (25) becomes 


(yay + Byr + 2p)0® — {p(B + lay + Bly + lyr + (y + B)zr}O 
+By=0 . (25a) 
If q = 0, equation (25) becomes 


6? — {(8 + l)ap + (y + lyr + (y + B)zr}® + Bre + yyr 
+ Byzp = . (250) 


We thus have a simple method of calculating the minimum reflux ratio. 
Equation (25) is used to find 6, and this value of 6 is then substituted in 
equation (22) to find R, which is the required minimum reflux ratio. 
Alternatively, S, the minimum reboil ratio, can be found from equation (23). 
Where y and z are the key components, 6 = }, = ¢, and lies between 1 
and 8. Where x and y are the key components, 0 = }, = ¢, and lies 
between Bandy. In both cases the value of 6 required from equation (25) 
is the one which lies between the relative volatilities of the key components. 

The method of calculation is illustrated by the following examples taken 
from the paper by Colburn.® 

Example 4. Third component lighter than the keys. 


paper Yp = 0-2; Zp = 0-2; q = l. 
» = 0-75; Yp = 0-25; ZD = 0; ty > 0; Yw = 9; Ze = I, 
pea 6B = 2. 


Equation (25) or (25a) becomes 


0-4 0-2 


— -+ — = 0 or 30? — 10-406 + 8=0 


1—¢ 
The desired value of 6 is that which is between 1 and 2. The solution is 


found to be 1-152. 
Equation (22) then gives 


= R+- land R = 0-643. 
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Example 5. Third component heavier than the keys. 
Lp = 0-2; yp = 0-2; zp = 0-6; gq = 0. 
Zp=—1; yp = 90; 2=p = 0; te =0; Ye = 0-25; zy = 0-75. 
y=10; 8 =2. 
Equation (25) or (25b) becomes 


2 0-4 | 06 
- _ il ae 
oe te + p= 1 OF Ot — 100 + 14-4 = 0. 


The desired value of 6 lies between 2 and 10, and is found to be 8-255. 
Then from equation (22), 


10 
10 — 8-255 


In a paper by, Mayfield and May * a method was given for calculating 
minimum reflux ratio which was based on the hypothesis that the minimum 
reflux ratio of a ternary mixture was the same as the minimum reflux ratio 
of two binary mixtures into which the ternary mixture was resolved. The 
hypothesis was not proved, but was demonstrated to give results correspond. 
ing to the results obtained by other methods for a number of examples in 
which the condition of the feed was g = 1 org = 0. It can, however, be 
proved correct for the general case by means of the equations which have 
been derived in this paper. 

Consider the case where the third component is heavier than the key 
components, which are x and y. The feed composition is xp, yp, zp. Let 
the ternary mixture be divided into two binary mixtures. The first binary 
mixture contains part of component 2, say zy’, and all of component y. 
The second binary mixture contains the rest of component x, say xp’ and 
all of component z. 

The feed composition for the first binary mixture is zp’ YP 





= R+1land R = 4-73. 





tp’ + yp’ Xp’ + yr 
The relative volatility of the two components is :. Applying equation (25) 


to this binary mixture, the term m. 4 corresponding to the third com- 





1 —6 
ponent, disappears. The equation becomes 
Y ° Xp’ 
Ry WP _ (1 — ghey’ + yr) (26) 
Y_ 1— 0 af - sll ita ileal 7 
6 


For this binary mixture, with x)= 1, equation (22) becomes 
fA 


= R’ + 1, and substituting this value of R’ in equation (26) gives 


Y 


Byp(R’ + 1) 
8B —(y — 8)R’ 


2x’ p(R’ ot 1) a 


= (1 —g)(z'y +yr)- - (27) 








ent © © = 





'S 
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Similarly for the other binary mixture, with relative volatility y, there is 
obtained the equation 


ap"(R" + +1) + op = (l= aer" +2)» 28) 


Now, by the hypothesis, if the minimum reflux ratios for the two binary 
mixtures are made equal to each other, they should also be equal to the 
minimum reflux ratio for the ternary mixture. 
Putting R’ = R” = R in equations (27) and (28) and adding them, we 
have, since 2’p + 2p» = Xp and tp + yp + zp = 1, 
Byr(R + 1) zp(R + 1) . 
ap(R +1 racine ie eee hme io: 8 
AR +1) + op tote q + (29) 
For the ternary mixture, the minimum reflux ratio is given by equations 
25) and (22). From the latter equation we have, since x» = I, 





yY _ a 
reer ty =R+loré= Rit 


With this value of 6, equation (25) becomes 


ty(R+1) , 2(R+1) _y_, 
—(y—8)R '1—(y—UR 


which is the same as per (29) derived from the two binary mixtures. 

A similar proof can be given for the case where the third component is 
lighter than the key components. For this case it is more convenient to use 
equation (23) with equation (25) to give the final result in terms of the 
minimum reboil ratio. 

It will be noted that equation (29) can be used to give a direct solution 
for R, instead of following the procedure previously described of first 
finding 6 from equation (25) and then finding R from equation (22). This 
latter procedure is the more convenient one, as it readily permits of the 
correct value of 6, and therefore of R, being selected. The use of equation 
(29) involves finding the various values of R which satisfy it and then 
deciding which value correctly represents minimum reflux conditions. 


xp(R + 1) +3 
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FRACTIONAL DISTILLATION OF MULTI-COMPON. 
ENT MIXTURES—CALCULATION OF MINIMUM 
REFLUX RATIO. 


By A. J. V. UNDERWoop, D.Sc., M.I.Chem.E., A.M.I.Mech.E., 
F.R.L.C., F.Inst.F. (Fellow). 


SuMMARY. 


By extending a method which has been applied to ternary mixtures, 
equations are derived by which the minimum reflux ratio for multi-component 
mixtures can be readily calculated when the fractionation between the key 
components is a sharp one, relative volatilities and molal reflux being 
assumed constant. Numerical examples are given to illustrate the method 
of calculation. 


THE calculation of minimum reflux ratio is an important problem. 
Several methods of calculation have been put forward such as those of 
Brown and Martin,! Brown and Holcomb,? Gilliland,? Hogan,‘ and (Col- 
burn.5 Some of those methods involve a considerable amount of calcula- 
tion and, as has been pointed out by Colburn,5 are open to criticism in regard 
to the basic assumptions made. The most convenient method for most 
purposes is that of Colburn,® but it admittedly rests on an empirical basis. 

By extending to multi-component mixtures the methods previously 
described by the author®’.in connection with ternary mixtures, the 
calculation of minimum reflux ratio can be made fairly simply for cases of 
ideal mixtures with constant relative volatilities and sharp separation 
between the key components. Where the relative volatilities vary in 
the column the same method of calculation can be used to give the minimum 
reflux ratio approximately. 

By way of example a mixture of four components (w, x, Y, 2) will be con- 
sidered, the relative volatilities, referred to component (z, being 8, y, $ 
respectively. It has been shown previously ® that an equation can be 
derived connecting the composition on any plate (wp, 79, Yo, Z) with the 
composition (w,, X,, ¥,, z,) on the next plate below it, namely 


buy, Y*1 BY, a: “ty 
Rate gt+y—eta—eti-g 
m(dw, + yt + BY, + Bi) 





(1) 


where ¢ is given by the equation 
as by cB 
oT y-$ 8-3 


As previously, for a rectifying column, m = 


“+ 








i— 
ee ee 
R+1 Frr R+1 


etc., and for a stripping column, m = oie whee 3 = RP + . —W 





~~ 


anda = — zy? 


tw 
wy? etc. 





/N- 
JM 


i? 


, 
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Equation (2) gives four values of ¢, which, in ascending order of magni- 
tude, are denoted by ¢,, $9, $3, $4. These values of ¢ give four equations 
when substituted in equation (1). Applying these to n successive plates 
and dividing any one equation by any one of the others, there are obtained 
a set of equations such as 











bw Yo Yo 20 
8— “eer ee me Rae cr 
dw = Yo PY _ + _* 
8—¢ y—b, B—¢d¢, 1— ge 

3w, 1 Yn BYn zn 

_ (#1) Fh “y= 8 * 4 IS * 
Oo me |, Me 
3—¢ y—¢, B-—¢ 1— dy 


Any two of the four values of ¢ will give an equation of this type, so that 
six equations are obtained in all. Denoting, for convenience, the expression 


5wg _Y%o , Yo | . E ‘ 
sh ty—b tag Tigo - «Om 


the six equations, of which equation (3) is the first, are 


E(x», $1) = (#1) E(x, $1) . E(x», $1) _ (#1) E(xp, $;) 

E(x, $2) Po! — E (Xn, bg)’ E(xq, $3) p;/  E(xp, 5) 
Hew $1) _ (#1) Blew $s), Elt~ $y) _ ($2) Sem $y 4) 
E(x, $4) dy — E(xy, $4) : E(x, $3) $3’ — E(xy, $3) 
E(x $2) __ ( y E(w $2). BZ» $s) _ ($s) E(xn $3) 
E(xo, $4) ps " E (Xn, $4) F E(x, $4) dy’ ~ Elan, $4) 
These six equations represent only three independent relations. The 

last three can obviously be derived from the first three. The three inde- 


pendent relations, combined with the equation 
w+r+y+z=1 
suffice to solve for the four variables, w, x, y, z. 

Under certain conditions the composition on a plate may be exactly 
equal to the composition on the plate above or below it, and a plate-to- 
plate calculation leads to an infinite number of plates. For such conditions 
of “‘ constant composition ’’ or “ limiting composition,” denoting such a 
composition by (g, h, &, 1), equation (1) gives 


g=m(’got+yh+pk+ . ..... (6) 
As shown previously ® 


en lee ees eee ee er eee | 
iat iy eke cy ee Hr ae or ee 




















There are four values of ¢, and consequently four limiting compositions 
which may be denoted by the suffix attached to the corresponding value 
of d. Thus g,, hy, k,, 1, are obtained by using ¢, in equations (6). 
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The expression H(z, ¢,) can be shown to become zero if gg, hg, ke, I, are 
substituted for w,, XZ», Yn; Zn. The condition for this is that 


























892 hg Bk ly =0 
3— ¢, nr 3 +3 — $y sf — $4 
or, from equations (6), 
as ‘ by ae cB 
(8 —#i)(8 — bo)  (y — bly — $2) — (8 — $108 — $2) 
d " 
Taam)" 
From equation (2) 
as by cB d 
5 — 4, eee 4 +; — $1 *% — 4 
A. 8 ee 
“S| Fk Foe tHe ; 
or 
as 4 by 4 cB 
(3 — $:)(8 — $2) © (y — oy — $2) (8 — $118 — $s) 


d 


ae ger = po» ale 





which is the same as equation (7). 

Similarly it can be shown that E(x,, ¢,) also becomes zero for the values 
93» hs, ky, 1, and gy, hy, ky, 1. It does not become zero for the values 
93, hy, k,, 1,, as may be seen by putting ¢, = ¢, in equation (7). All the 
terms on the left-hand side are positive for a rectifying column and negative 
for a stripping column, and the sum of them can in neither case be zero. 

Similar relations can be obtained for the other expressions similar to 
E(x,, ¢;). The various equations are satisfied as follows. (The suffix n is 
here dropped for the sake of generality.) 


E(x, $1) = 0 by (go, he, ke, Ip), (9g, hg, Hg, tg) and (gq, hy, ky, 1) 
E(x, ¢) =0 by (91, h,, ky, l,), (93; hs, ks, Is) and (Ja, hy, ky, 1.) 
E(x, $3) = 0 by (91, Ay, &,, 4), (92, hg, ke, 1.) and (94, hy, ky, 1) 
E(x, $4) = 0 by (93, 4, ky, 11), (Ga> he, ky, I) and (gs, hg, ks, 13) 


If any three of these equations are simultaneously satisfied, they have one 
common solution. 


Thus E(x, $4) was E(x, $3) _ E(x, $4) =0 . 


have (9g, hg, ky, 1) as the common solution. The suffix attached to g, h, k, | 
is the one which does not appear in the values of ¢ involved in the three 
equations. 

If any two of the equations are simultaneously satisfied, they have two 
common solutions. 


Thus E(zx, $3) = E(x, $,) = 0 





W 


noc 
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have (93, 4y, &y, 1) and (gq, hy, ky, 12) as common solutions. Here again the 
suffixes attached to g, h, k, 1 are those which do not appear in the values 
of ¢. 

fie significance of the four limiting compositions can be shown as follows, 
a rectifying column being considered by way of example. It is readily 
shown, as in the case of ternary mixtures, that putting wy = wp, % = Zp, 
etc., makes any of the expressions on the left-hand side of equations (4) 
equal to (1). We then have, by selecting the three appropriate equations, 


E(%n, $2) _ ($1)". E(u $3) _ ($1)". E(t $4) _ ($1)" 
E(2p, $)) (f) ? E(x; $,) () , E(x,, $1) ($!) 
Now ¢,<¢_.<¢3<¢, and therefore, when n becomes infinite, 


E(x, $2) = E(xp, $3) - E(x, $4) = 0. F: 

These equations are satisfied by (g,, 4,, &,, /,), which, in addition, does not 
make the denominator E(zx,, ¢,) zero. Thus a calculation proceeding down 
the column from the top finally reaches the limiting composition 
a» Ay Ay 4). 

If we put n = — N, where N is positive, the composition can be obtained 
for an imaginary plate located N plates above the top of the column. From 
equations (4) we have 


Beam dh) (fi)", Ben dd) _ (fs)", Ben de) _ (fs)! 














$s ; E(z_y, $4) d,s 
When n becomes infinite 
E(x_y, ¢;) - E(x_y, $2) = E(x_y, $3) = @ 


and the solution is (g,, h4, k, |,), which is the limiting composition reached 
when a calculation is carried upwards from the top of the column. 

The other two limiting compositions (g,, hg, ky, 1.) and (g5, hg, ks, ly) have 
particular practical significance when certain components are only present 
in small amount in the top product. 

There are three cases to be considered :— 


Case I. zis the heavy key component and is present in the top 
product in small amount. 

Case II. yis the heavy key component and both y and z are present 
in the top product in small amounts. 

Case III. x is the heavy key component and z, y and z are present 
in the top product in small amounts. 


For case I, since d in equation (2) is approximately zero, ¢, is approxi- 
mately equalto 1. Substituting this value in equations (6) gives g,,h,, k,l, 
the last of these by difference. These are the values of the limiting com- 
position reached by calculating downwards from the top of the column with 
all four components in the calculation. 

For case II, since both d and c in equation (2) are approximately zero, 

, = 1 and ¢, = 6 approximately. Using the value ¢, = 1, the values 
of g,, Ay, ky, 1, can be calculated when the values of d and c, although very 
small, are specified. The method of calculation is similar to that given for 
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Example 3 of the previous paper.” Using the value ¢, = 8 in equations (6) 
ives go, hg, ka, 14. 1, has a very small negative value which, as explained in 
connection with ternary mixtures,’ may be neglected for practical purposes, 
k, is indeterminate in equation (6) when c = 0 and ¢ = 8, but can be found 
by difference, or by the method of approximation described in connection 
with ternary mixtures. This limiting composition thus contains the three 
components w, x, y, and is the one reached by calculating down from the 
top of the column using these three components. It can also be calculated 
by the usual method based on the small amount of the heavy key component, 

For Case III, since d, c, and 6 are all small in equation (2), d, =1, 
¢, = 8, and ¢,= y approximately. Using ¢, = y in equation (6) gives 
kg =l,=0. After finding g,, then h, is found by difference. The limiting 
composition contains only the two components w and zx, and is the one 
obtained by calculating down the column from the top using these two 
components. It represents the limiting composition that would be reached 
with a binary mixture. 

Thus, when making a calculation downwards from the top Op hy, k,, l) 
is reached if z is the heavy key. If y is the heavy key, (go, hg, kg, I.) is 
reached. If x is the heavy key, (93, Ag, kg, ly) is reached. 

For a rectifying column, superimposed on a stripping column, minimum 
reflux conditions prevail (in the rectifying column) when a calculation made 
downwards from the top and a calculation made upwards from the feed- 
plate both meet at the appropriate limiting composition which, as shown 
above, depends on which are the key components. It is possible to specify 
what composition or range of compositions can exist at the feed-plate if the 
downward and upward calculations are to meet in this way. 

If z is the heavy key, the composition at the feed-plate will be 
(93, 4, ky, 1,), and the downward and upward calculations actually meet at 
the feed-plate. The composition at the feed-plate is thus defined by the 
three equations 


E(x, $2) = 90; E(x, ¢3)= 0; E(x,d)=0 . . . (8) 
If y is the heavy key, the limiting composition (gg, hg, ky, 1.) is reached by 


the downward calculation. If the upward calculation is to reach this. 


point, H(x,, ¢,) must not be equal to zero. (Here wp, Zp, Yn, Z is taken to 
be the feed plate composition.) From equation (1) it will be seen that, if 
E(x, $,) is zero for any given plate, it will also be zero for all successive 
plates. From equations (4) we have, for a plate 0 which is n plates above 


the feed plate, 
E(x, $1) _ (fy E (2p, $1) 


E(x, $2) $e " E (xn, $2)" 


When n becomes infinite, H(x9, 4,) becomes zero. At (gg, hg, kg, Jy), in 
addition to E(x», ¢;) = 0 we must also have E(ao, $3) = Oand E(x», ¢4) = 0. 
If these two equations are satisfied by the feed-plate composition, they will 
also be satisfied by the compositions on all plates above it, as will be seen 
from equation (1). Thus, where y is the heavy key, there is a possible range 
of compositions at the feed-plate, and this range is defined by the equations 


(9) 











~ E(x, $5) = 90; E(x, d) =90 
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From these equations it is also clear that the upward calculation will not 
pass through (93, hg, kg, Uy) Or (Gq, Ags Ky, U4). 

If x is the —. key, the limiting composition reached by the downward 
calculation is (93, hg, ks, 1s). This composition must also be reached by the 
upward calculation. Then E(x», ¢3) must not be zero. 

From equations (4), 


ace ¢}) ai E(x, $;) d E(x, $2) $s E(zx,, 2) 
E(x, $3) = (fy. " E(an, $5) vet E(x, $3) =( ) " E (2p, $3) 


When n becomes infinite, H(29, ¢;) = Oand E(x», d,) = 0. At (gg, hg, ky, Is), 
we must also have E(x, ¢,) = 0. If the feed-plate composition satisfies 
this equation, the upward calculation will reach (95, hy, ks, 1). Thus, where 
x is the heavy key, the possible range of compositions at the feed plate is 
defined by the equation 
E(z,¢,)=0 . . . o ¥- are 

From equations (8), (9), and (10) it will be seen that the composition at 
the feed-plate, for minimum reflux conditions, is defined by three, two or 
one equations according as z, y, or x is the heavy key. 














Conditions for minimum reflux in the stripping column can be defined in 
asimilar manner. When calculating upwards from the reboiler, the limiting 
composition (gq, ha, kq, ly) is reached when w is the light key. When z is 
the light key, (g5, hg, ks, ls) isreached. When yis the light key, (Ja he» ke, 1) 
is reached. By calculating downwards below the reboiler. (93, hy» ky, le) 


is reached. 
For the stripping column it is convenient to write » instead of ¢. Also 


Yy<bg<by<ty. Let (wo, x9, Yo 2) refer to the feed-plate and (wy, Xn, Yn» Zn) 
to a plate n plates below the feed-plate. When y is the light key, 
(ga, has Kg: lg) is reached by calculating upwards from the bottom, and must 
also be reached by calculating downwards from the feed-plate. Then 
E(Zp, ¥2) is not zero. 
From equations (4) 
E(xn, Ys) Ls (%2 , E(x, Ys) 1 E(xn, V4) _ = (4)". E(x, Ys) 
E(x; Ye) 3’ " E(x; be) E(Zn; be) YW E(x», be) 
When n becomes infinite, E(z,, }3) = Oand E(z,, ¥4) = 0. For (go, he; ke, le) 
to be reached it is also necessary that E(x,, ¥,) = 0, and this will be the 
case if the feed-plate composition satisfies the relation 
E(z,¥,)=0 .. .  * a a 
If x is the light key, (93, hgs kg, lg) has to be senile by calculating down 
from the feed-plate. 
E(x, 3) is not zero. From equations (4) 
E(Zn, Yq) _ (3) E(xo, Ya) 
E(x», 3) E(x, bs)" 











be 
When n becomes infinite, E(z,, ¥,) = 0. In addition, for (g5, hg, kg, Is) it is 
necessary that E(x,, ¥,) = Oand E(z,, v2) = 0. This will be the case if the 
feed-plate composition satisfies the relations 


E(x, $3) = Oand H(iz,¥)=0 . . . . . (12) 


ae 
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If w is the light key, the limiting composition (g,, hy, k,, ly) is reached at 
the feed-plate and the relations to be satisfied are 


E(x, ¥,) = 90; E(x, a) = 90; E(x, 4.) = 0 . . ~ (13) 


Relations have now been derived which specify the composition or range 
of compositions at the feed-plate if minimum reflux conditions prevail in 
the rectifying column or in the stripping column. If minimum reflux 
conditions prevail throughout the whole column, the conditions for feed. 
plate composition derived from the rectifying column and those derived 
from the stripping column must be satisfied simultaneously. Collecting 
together the results which have been obtained, the conditions to be satisfied 
for the various separations are shown in Table I. 








TasBLeE I. 
= = y | mom rectifying column. From stripping column. 
y z {Pl $2), $¢) = way d ae. ¢s) = E(x, $,) = 0 , 
2 y E(z, re = 0; E(x, ¢,) = 90 E(x, $,) = 0; E(x, ¥,) = 0 
w x E(x, $,) = 0 (Fe ks, p Bre a i "7 

















In each case there are four equations to be satisfied. Each of these 
equations is of the type 


dw yx By z 
=0 
or ae 1 Me em Sk ee * 


and involves only three independent variables such as = , and 2. (In 











addition, there is, of course, the equation w + 2 + y + z= 1 for finding 
the actual values of the components from their ratios.) There must, 
therefore, be a relation between the four equations which reduces them to 
three independent equations. A relation which obviously meets this 
requirement is that any two of the values of ¢ or | involved in the equations 
should be equal. 

From equation (2) it is readily shown that, for a rectifying column where 
a, b, c, d are positive 


0<¢,<1; 1<¢.<8; B<¢3<y; y<¢y<3. 
Similarly for a stripping column where a, b, c, d are negative 
1<$,<B; B<de<y; y<dg<8; 8<yy. » 


Thus none of the values of ¢ can be equal and none of the values of } can 
be equal. To satisfy the conditions for minimum reflux it is therefore 
necessary that one of the values of ¢ involved should be equal to one of the 
values of | involved. 

For the case where y and z are the light and heavy keys respectively, 
}, must then be equal to ¢, or ¢, or ¢,. From the ranges of the values of 
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¢ and y, it is obvious that the only possibility is ¢. = 4,. Both these lie 
between 1 and 6—that is, between the relative volatilities of the key 
components y and z. The limiting composition in the rectifying column in 
this case is (g,, hy, k,, 1,) which satisfies E(x, $,) = 0. The limiting com- 
position for the stripping column is (g», he, ke, 12) which satisfies H(x, »,) = 0. 
Since ¢2 = $, both limiting compositions satisfy the same equation. 

Where « and y are the light and heavy keys, 4, or 4, must be equal to 
$, or dy. The only possibility is ¢, = Y, as both lie between 6 and y, the 
relative volatilities of the key components. Here again it can be shown, 
as before, that the limiting compositions for the rectifying and stripping 
columns both satisfy the same equation E(x, ¢,) = 0 or E(x, },) = 0. 

Where w and z are the light and heavy keys, ¢, must be equal to 4, or 
J, or ¥g- Thus ¢, must be equal to 4, as both lie between y and 8. 

The condition for minimum reflux in both the rectifying column and the 
stripping column is therefore 


ee a 


where n and (n + 1) represent the numbers to be attached to the values 
of ¢ and » when counting up from the lowest values, and also represent the 
numbers to be attached to the heavy and light keys components when 
counting up from the least volatile component. Also ¢,,, and , have a 
value which lies between the relative volatilities of the key components. 

When equation (14) is satisfied, a stepwise calculation starting from the 
reboiler will reach a limiting composition in the stripping column and an 
infinite number of plates will be required to pass through this composition. 
The calculation will then pass up through the feed plate and will reach a 
limiting composition in the rectifying column and an infinite number of 
plates will again be required to pass through this composition to reach the 
composition at the top of the column. 

Equation (2) for the rectifying column and the corresponding equation 
for the stripping column thus have a common root, which will be denoted 
by 6. ‘This root has a value lying between the relative volatilities of the key 
components. Substituting for a, 6, c, d in equation (2), the appropriate 
values for a rectifying column and for a stripping column, there are obtained 
the two equations which define 6 





8wp yXp BY Zp _ 

so ty-0 s—etr—enm Ft! dtd sp. Seca 
and 

suw | yt byw | wm __g 

s—ot+y—otp—oti—o- S tin 0 (16) 
Now Pw, + Www = Fwy; Pry + Wry = Fry; ete. 


where wy, <p, etc., denote the composition of the feed. 

Also SW = RP + qF — W =(R+1)P — (1 — q)F since W = F — P. 
Then multiplying equation (15) by P and equation (16) by W, and adding, 
we obtain 








dwp yur SYP  —— oe 
a fat baal Si a be 
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When the feed is liquid at boiling-point, g = 1. When the feed is vapour 
at dew point, g = 0. 

From equation (17) is found the appropriate value of @ which is the one 
lying between the relative volatilities of the key components. On sub. 
stituting this value of 6 in equation (15), the minimum reflux ratio js 
obtained. Similarly the minimum reboil ratio can be obtained from 
equation (16). 

When R and S have thus been found, equations (19) and (20) can be 
used to find the other values ¢ and 4 required in the equations for the feed. 
plate composition listed in Table I. These equations can then be used to 
solve for the feed-plate composition. 

In equations (15), (16), (17) the relative volatilities are all referred to the 
least volatilecomponent. The relative volatilities can, however, be referred 
to any component. Thus, if y is the heavy key, and it is desired to refer 
the relative volatilities to it, equation (17) takes the form 








=Il1—-q ° ° . (18) 


Equations (15) and (16) can be written in the same form. The variable is 


now —, which may be written 6’. Its value lies between 1 and rie, 


between the relative volatilities of the key components. The value obtained 
for the minimum reflux ratio is obviously the same whichever form of the 
equations is used. 

The general equations corresponding to equations (15) and (16) for 
conditions other than minimum reflux are 








2D vs 
Poet Rti. e - 9) 


(20) 





By differentiating these equations it will be seen that ¢ increases as R 
increases, while ) decreases as S increases. FR and S must obviously increase 
together. When R and S become infinite, the roots of both equations (19) 
and (20) are obviously 1, 8, y, 8. As R decreases from infinity to a finite 
value, the values of ¢ decrease and the values of 4 increase. As R is 
gradually decreased there will come a point when a value of ¢ and a value 
of » (lying between the relative volatilities of the key.components) will 
become equal, and this will correspond to minimum reflux conditions. 

The analysis has been given in detail for a four-component mixture. It 
can be made in exactly the same way for a mixture of any number of com- 
ponents, and equations exactly similar to equations (15), (16), (17) are then 
obtained. For instance, for a six-component mixture there will be six 
equations to be satisfied at the feed-plate for minimum reflux conditions. 
The ratios of the components constitute five independent variables. Two 
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of the equations must be the same, and this condition is found to be 
equation (14) 
pn ce = Vn 


where n and (n + 1) refer to the heavy and light key components and the 
common value of ¢ and y lies between the relative volatilities of those 
components. The relation is quite a general one, and also holds good for 
binary mixtures. 

The root of equation (17) which is normally required is the one lying 
between the relative volatilities of the key components. In some cases it is 
useful to obtain the other roots. For instance, if a mixture of four com- 
ponents is to be separated into the four pure substances, three fractionating 
columns are required. There are five different ways of using these three 
columns to produce the four pure substances, as illustrated diagram- 
matically by Thormann.® In such a case it is useful to calculate the heat 
requirements for the different ways of effecting the total separation. 
Equation (17) is then used to find the three values of 6 corresponding to the 
cases where one, two, or three components are taken as the top product of 
the first column. Equation (15) is then used to obtain the corresponding 
minimum reflux ratios, it being noted that equation (15) varies with the 
composition of the top product, while equation (17) does not. The mini- 
mum reflux ratios for the other two columns, which may be separating 
either ternary mixtures or binary mixtures, are calculated similarly for the 
various cases. 

The relative ease with which the calculation of minimum reflux ratio for 
multi-component mixtures can be made is illustrated by the following 
numerical examples. 


Example 1 (from Colburn 5) 
wand x are the key components 


Wp = Ly = Yp = ty = O25; wp = | 
8=8; y=4; B=2. gq=1. 


Equation (17) gives ase 


> « 
Aas 


y 5 
i ee es. 2 ee 
§—6-' 4—6 '3-—e' 1-0 


The value of 6 required lies between 4 and 8. It is found by trial solution 
to be 5-58. 
From equation (15), 
2 
8 — 5:58 
Example 2 (from Colburn 5) 
As in Example 1, but with x and y as key components 


=R+1 «R=231 


Wp = Lp = 05. 


Equation (21) is used again, the value of ® now required being that 
between 2 and 4. It is found to be 2-556. 
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From equation (15), 


4 4 2 
8 — 2-556 © 4 — 2-556 





=R+1 .«. R=1-12 


Example 3 (From Colburn 5). . 
As in Example I, but with y and z as key components 


Wp = Xp = Yp = 0-3333. 


Equation (21) is used again, the value of @ now required being that 
between 1 and 2. It is found to be 1-196. 
From equation (15) 
2-667 1-333 4 0-333 
8 — 1-196 ' 4—1-196 * 2—1-196 


Example 4 (from Gilliland %), 








Relative 
Distillate. Residue. volatility 
to o-cresol. 








Phenol . . : ‘ “3: QE 0-0524 
o-Cresol . ‘ . ‘ oD D 0-199 
m-Cresol . ‘ 3 . 0-449 
Xylenols : : ; , 0-224 
Residue . : 4 ‘ ; 0-075 

















Here there are five components. Relative volatilities are referred to 
o-cresol. Instead of recalculating them relative to the heaviest component, 
an equation of the type of equation (18) can be used with the variable 
changed to 6’. An appropriate term for the fifth component is brought into 
the equation. (gq = 1) 


1:26 x 0:35 0-15 | 0-663 x 03 0-304 x O-15 
126-0 ‘i-#t 063—0 ' 0304 —0 
0-087 x 0-05 _ 9 
0-087 — 0’ 


The value of 6’ required lies between 1 and 1-26. By trial solution it is 
found to be 1-0798. Successive approximations and interpolations are not 
laborious, as, when an approximate solution has been obtained, only the 
first two terms in the equation change appreciably with small variations in 
the assumed value of 6’. 
The equation corresponding to equation (15) for the modified variable 
0’ is . 
1-26 x 0-95 1 x 0-05 
a — ye * 1-2 aoe. 


Gilliland * states that “‘ detailed stepwise calculations indicate that the 
true minimum reflux ratio is approximately 5-2.” The difference between 


this and the figure calculated above is thus 3-4 percent. In this example 
the amounts of the heavy key in the distillate (0-05) and of the light key in 











FRACTIONAL DISTILLATION OF MULTI-COMPONENT MIXTURES. 625 


the residue (0-0524), though small, are not negligible. This may affect the 
accuracy of the method used in this paper, which is based on negligible 
amounts of the heavy and light keys in thedistillate and residue respectively. 
On the other hand, plate-to-plate calculations are not easy to carry out 
with a high degree of accuracy, owing to the possibility of small errors being 
cumulative. 

This example may be used to illustrate the effect on the minimum reflux 
ratio of a comparatively small change in the relative volatilities of the key 
components when these are fairly close to each other. Assume that the 
relative volatility of phenol to o-cresol is 1-28 instead of 1-26, all other data 
remaining unchanged. Equation (22) then gives 6’ = 1-0856 instead of 
10798, and R is found to be 4-67 instead of 5-02. Thus a difference of 
about 14 per cent in the relative volatility makes a difference of about 7 per 
cent in the minimum reflux ratio. This is approximately the same differ- 
ence as for a binary mixture where R = —;, : Dz and the effect on R is 

Gk F 
proportional to (« — 1). It appears that a high degree of accuracy in the 
method of determining minimum reflux ratio is unwarranted unless the 
relative volatilities, particularly of the key components, are known to a 
correspondingly high degree of accuracy. Especially is this the case when 
the relative volatilities of the key components are not greatly different. 


Example 5 (From Jenny *). 





Relative volatilities at— 
Distillate. 





63° F. 175° F. 300° F. 





0-434 _ 514-5 
0-150 — 100-3 
0-411 0-010 34-1 
0-005 0-417 10-69 
_ 0-274 3°35 
— 0-299 l 























The feed is 66 per cent vapour and 34 per cent liquid, so that g = 0-34. The 
temperatures at the top and bottom of the column are 63° F. and 300° F. 
respectively. 'The temperature of the feed is 175° F. 

Jenny gives relative volatility data for all components at the feed tem- 
perature, but only for some components at the top and bottom tempera- 
tures. The additional data used have been taken from those published 
by Kirkbride.° The relative volatilities vary considerably through the 
column. For the key components the variation is from 1-86 at the bottom 
to 3-19 at the top. 

The method of calculation described in this paper is based on constant 
relative volatilities. To apply it to this example it appeared a reasonable 
approximation to take relative volatilities at a temperature midway be- 
tween the top and bottom temperatures. Actually the feed temperature 
(175° F.) has been taken instead of the true mean (181-5° F.), owing to the 
data being readily available for the former temperature. 
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Using these data in equation (17) with two additional terms, since there 
are six components, we have 


100 x 0:26 24-6 x 0-09 4 10 x 0-25 4 4-85 x 0-17 
100 — 6 24-6 — 0 10 — 0 4-85 — 0 
2-08 x 0-11 0-12 


ws Pnselacaagyl Dec 2: 
2-08 — 6 be es oF .. a 





The value of @ required lies between 4-85 and 10. It is found by trial 
solution to be 6-73. Substituting in the equation corresponding to equation 
(15) gives 

100 x 0-434 | 246 x 0-150 10 x 0-411 _ R+1 

10 —673 * 46—673 ‘10-673 "> 
from which R = 0-93. 

Jenny ® obtained a value of 0-95 and Colburn a value of 0-96. The 
agreement with these values is quite good. The general validity of the 
basis assumed in the above calculation—namely, taking relative volatilities 
at the mean temperature in the column—requires, however, to be checked 
by a larger number of examples. 

The limits between which the true minimum reflux must lie can be found 
by making two calculations, one based on the relative volatilities at the 
top of the column and the other based on the relative volatilities at the 
bottom. The first calculation assumes that the relative volatilities through- 
out the column are equal to the high values at the top and obviously gives 
too low a value for the minimum reflux ratio. Similarly, the other calcula- 
tion obviously gives too high a value. 

The value of R found by the first calculation is 0-53 and that found by the 
second calculation is 1-15. The mean of these is 0-84, compared with the 
true value of 0-95 or 0-96. As this mean value is obtained by averaging two 
substantially different values, the approximation can only be regarded as 
a rough one. 

For many practical purposes the higher limit for the minimum reflux ratio, 
obtained by taking the low relative volatilities at the bottom of the column, 
provides a useful figure for guidance and, as shown, it can be quite readily 
calculated. 
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THE TENDENCY TO SMOKE OF ORGANIC 
SUBSTANCES ON BURNING. PART I. 


By A. E. Cuarke, Ph.D., T. G. Hunter, D.Sc., A.R.T.C., M.I.Chem.E. 
(Fellow), and F. H. Garner, O.B.E., Ph.D., F.R.1.C. (Fellow). 


SuMMARY. 


The British incendiary bomb, used in very large quantities by the R.A.F. 
for the bombing of German towns and cities, was essentially a 30-lb. bomb 
filled with a special benzol gel together with white phosphorus. The benzol 
gel filling produced, on burning, a large amount of black carbon smoke which 
obscured the target, and resulted also in an appreciable portion of the filling 
being wasted as unburned carbon. In addition, copious white smoke pro- 
duced by the burning phosphorus increased the obscuring effect over the 
target. The possibility of replacing this benzol gel-phosphorus combination 
with a smokeless filling of a satisfactory nature was therefore investigated. 

Towards this end the smoking tendency of a large number of organic 
compounds was assessed by flame-height measurements in a special lamp 
based on the I.P. smoke lamp. 

A burning organic substance has a flame-height at, and above which, 
smoking occurs, and this height is a measure of the tendency to smoke. A 
new form of lamp was devised to measure flame-heights, from about 9 to 
450 mm, of liquid compounds burning freely in air. 

A wide range of hydrocarbons, alcohols, ketones, esters, and nitro-com- 
pounds was examined—115 compounds in all. In general, a compact 
molecule was found to give a smoky flame. The order for increasing 
tendency to smoke for hydrocarbons is: n-paraffins (in which increased 
chain length or chain branching gave increased smoke), naphthenes, olefines, 
and aromatics (in which appreciable aliphatic side chains on the benzene 
ring appeared to give €o marked reduction in smoke). 

In general, increased oxygen content of an organic compound resulted in 
decreased smoking tendency and compounds, such as methyl acetate, 
containing high percentages of oxygen only smoked at very large flame- 
heights. Some compounds, such as allyl alcohol, although having appreci- 
able oxygen contents, had relatively high smoking tendencies, due to the 
nature of the carbon-hydrogen portion of the compound. 

Of the aliphatic alcohols, the tertiary compounds were more smoky than 
the primary compounds. This also applied to nitro-paraffins. For each 
set of isomeric aliphatic esters, the flame-height at which smoking began 
increased with the chain length attached directly to the carboxylic carbon 
atom. 

At equal oxygen content, the general order for increasing tendency to 
smoke was: n-primary alcohols, n-primary nitro-paraffins, propionates, 
acetates, lactates, and formates, although the order varied slightly for 
different oxygen contents. 


WHEN an organic substance burns under fixed conditions there is a 
particular critical flame-height at and above which smoking occurs. It 
would appear that the tendency to smoke is determined by the amount 
of oxygen (required by the flame on the one hand and by the amount of 
oxygen) available to the flame on the other. At the critical flame-height 
the two quantities are equal. 

It has been observed by Bancroft ! that a luminous flame can be regarded 
as a colloidal suspension of carbon in a gaseous medium. Further, 
Minchin 2 suggests that through coagulation the charged carbon particles 
vary in size through the flame, being smallest towards the base and largest 
at the tip. Thus, in certain instances this coagulation may proceed to 
such an extent that the particles at the tip of the flame are too large to be 
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burned in the prevailing conditions of temperature and oxygen supply, 
and under these conditions smoking occurs. This coagulation rate, how. 
ever, depends on many factors, such as the concentration, the initial gize 
and the initial charge of the carbon particles, the rate of loss of charge, the 
mean flame temperature, and the composition of the gaseous medium, 
which, it appears, in their turn depend on the molecular structure of the 
compound. 

The amount of oxygen required by the flame from a compound under 
fixed conditions of burning and air supply will thus depend on : 


(1) The molecular structure, which determines the initial state and 
the properties of the carbon particles when the fuel is decomposed ; 

(2) The rate of burning of which flame volume is a function ; 

(3) The flame-height, which not only influences the final particle 
size of the carbon, but is also intimately connected with flame volume, 
and hence with rate of burning. 


On the other hand, the amount of oxygen available to the flame will 
depend on : 


(a) The oxygen content of the compound ; 

(b) The area of flame surface in contact with the air ; 

(c) The extent of dilution of oxygen at the flame surface by gases of 
combustion. 


Thus the flame-height at which smoking begins is intimately related to 
the tendency to smoke, and hence this critical flafhe-height can be used as 
a measure of the smoking tendency of the compound. Kewley and 
Jackson * have described a method of evaluating the tendency to smoke of 
kerosine, which is now employed as a standard method of assessing this 
property (I.P. Standard Method 57/45), in which the maximum height of 
flame that can be obtained just short of smoking is measured in milli- 
metres. The flame-height thus obtained is, in the standard method, used 
as a measure of the tendency to smoke. 

Minchin * employed with this flame-height lamp a slightly different 
method of expressing the value “ tendency to smoke ”’ using the formula : 


S, = 320/h 
where S, = smoking tendency 
h maximum flame-height in millimetres. 


Further, Terry and Field * have described an improved form of the Davis 
Factor lamp in which flame-heights of up to 102 mm could be measured. 


A New Smoxe Port Lamp. 


The maximum flame-height capable of measurement in the I.P. lamp is 
50 mm, but with the improved Davis Factor lamp the maximum measurable 
flame-height is 102 mm, while the smoking flame-height of a good kerosine, 
with this lamp, exceeds 75 mm.5 

A lamp was required for this work in which the flame-height could be 
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varied and measured over a wide range without the use of a wick. 
Preliminary experiments indicated that, when a substance burned, an 
increase or decrease in burning area resulted in a corresponding increase or 
decrease in flame-height. It was felt that this principle might be used to 
vary the flame-height in a lamp, and the method finally adopted for 
changing the burning area is illustrated in Fig. 1. A shallow conical 
metal vessel was connected by oil-resisting rubber tubing to a liquid 
reservoir as shown. When liquid was introduced, the level in the reservoir 
determined the level in the cone, and hence the exposed liquid area in the 
cone. ‘Thus, on burning the liquid in the cone, a variation in the height of 
the reservoir produced a corresponding variation in the burning area, and 
hence in the flame-height. 

It was apparent that for any accurate observation of flame-height a 
very steady flame was essential. With the large flames that could be 
produced when using the cone burner, however, it was found that even in 


RESERVOIR 
~3— POSITION 2 


AREA 2 


RESERVOIR 
POSITION | 


Fig. 1. 
THE PRINCIPLE OF THE CONE BURNER. 


the absence of draughts, bad flickering occurred. This difficulty was 
overcome by enclosing the flame and using the draught from a glass chimney 
to draw the air necessary for the combustion through two fine copper 
gauzes (120 mesh) before reaching the burner (see Fig. 2). Thus the air 
flow was stream-lined, and a steady flame resulted. 

As the determination of the exact point at which smoking began was 
somewhat difficult, it was decided to mark the position of the tip of the 
flame by a wire hook, which could be raised or lowered in the chimney by 
means of a screw. This wire passed through a hole in a copper gauze 
(18 mesh) at the top of the chimney, which damped any vibration set up 
when the wire was moved. The height of the hook was shown on a metre 
rule, the position of which was adjusted so that zero flame-height was 
indicated when the wire just touched the lowest point in the burner. 

Preliminary experiments with this lamp were carried out by burning 
both highly inflammable substances and those substances such as cresol 
which burned only with difficulty. With the latter type of compound it 
was found that the burner had to reach a high temperature before a steady 
flame was produced. Thus the burner was made from thin brass (specific 
heat 0-092 cal/gm), which gave it a low heat capacity. This alone, however, 
was unsatisfactory, since substances of high volatility tended to boil and 
produce an irregular flame. This difficulty was overcome by fitting to 
the base of the burner a cold-water-circulating jacket which could be used 
when burning these volatile substances. 
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A diagram showing constructional details of the lamp finally employed 
is given in Fig. 2, and a photograph is shown in Fig. 3. In Fig. 2 


A is the conical burner; 
B is the chimney ; 
C is the glass collar. 
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THE SMOKE-POINT LAMP. 





It is important that A should be of a larger diameter than B, and the 
distance between the base of B and the top of A is adjustable, so that flow 
is in the laminar regime, and thus is steady. With this lamp flame- 
heights were measured over the range of 9 mm to 450 mm. 
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Fic. 3. 
THE SMOKE-POINT LAMP. 
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Procedure. 


Having removed the chimney and glass collar, the reservoir, connection 
tube, and burner were washed with the liquid to be used and then filled to 
a suitable level. After igniting with a bunsen burner (some liquids 
required heating before they burned), the collar and chimney were replaced, 
and the liquid was allowed to burn for 3 minutes to obtain fairly steady 
conditions. 

The flame was then raised or lowered, by movement of the reservoir, to 
its smoke-point, which was indicated by the movable wire. It was found 
that each flame formed a long red “ tail ’’ just before smoking began, and 
in these tests the smoke point was taken as the height at which this “ tail ” 
just disappeared. A number of observations were made for each substance, 
and the mean value taken. The experimental error is considered to be not 
greater than +3 per cent. 

The smoke points of a wide variety of organic substances (115 in all) 
were determined in the lamp, and were used as a measure of the tendency 
to smoke—a high smoke point indicating a low tendency to smoke, and a 
low smoke point indicating a correspondingly high tendency to smoke. 

The compounds examined in the lamp were : 


(1) Hydrocarbons and hydrocarbon mixtures. 
(2) Oxygen containing compounds, including : 


(a) Alcohols. 

(b) Glycol and glycerol compounds. 
(c) Nitro compounds. 

(d) Ketones. 

(e) Esters. 

(f) Carboxylic acids. 


Results. 
(1) Pure Hydrocarbons. 


The smoke points of pure hydrocarbons as determined in the lamp are 
given in Table I. 

Of the three normal paraffins examined, n-pentane and n-heptane have 
similar flame-heights with cetane somewhat lower, showing a gradual 
increased tendency to smoke with increased chain length. Isopentane 
with a branched chain is more smoky than n-pentane. The flame-heights 
of n-hexane and n-octane must be close to those of n-pentane and n- 
heptane (i.e., roughly 15-16 cm), and 2 : 2-dimethylbutane with a highly 
branched chain is much more smoky, and similarly with the branched-chain 
octanes. Thus the branching of the carbon chain in a compound increases 
the tendency to smoke, and the results suggest that more than one branch 
produces a greater effect. The flame-heights of the three branched-chain 
octanes indicate that the positions of the carbon branches have little or no 
effect on the smoking properties. 

As the only olefines available were probably mixtures, no reliable con- 
clusion can be drawn, except that the introduction of the double bond, 
and probably some branching, make these compounds very smoky. 
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Benzene, having a high C/H ratio, is very smoky indeed, and the building 
up of this effect may be traced from n-hexane through cyclohexane to 
cyclohexene and benzene. A similar effect may be noted with decalin and 
tetralin. Results 14-17 for toluene, xylene, di-isopropyl benzene and 
p-cymene indicate that the introduction of aliphatic side-chains in the 
benzene nucleus has little effect on the smoky nature of the benzene flame, 

It is interesting to compare these results with those obtained by 
Minchin.* This investigator was able to show that the flame volume was 


Taste I. 
Hydrocarbons. 


Smoke point 
flame height, 
cm. 





. n-Pentane . 
. Isopentane : ‘ 
2 : 2-Dimethy]! butane 
n-Heptane . , . : : 
Iso-octane (2 : 2: 4-trimethyl pentane) | 
2:2: 3-Trimethyl pentane (84%) F 
2:3: 4-Trimethyl pentane (82%) 
Cetane . ‘ i . ° 
. Pentene 

. Hexene 

- Heptene . 

- Decene 

. Benzene . 

. Toluene , ; 

. “* Xylene ” : ‘ 

. m-Di-isopropyl benzene 

. p-Cymene . . 

. cycloPentane 

. cycloHexane 

. cycloHexene 

. Methylcyclohexane 

. Dicyclohexyl 
23. Pinene 

. Tetralin 

- Decalin 
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related to the molecular volumes of combusted products and of oxygen 
required for combustion by the formula : 


1 
prhgthe 


where V = flame volume; 
y = molecular volume of oxygen required for complete com- 
bustion ; 
x = molecular volume of combusted products ; 
K, and K, = constants. 
Since the flame, in the lamp he employed, was conical, then , = : x a 
constant, and hence hf (the smoke point flame-height), for any substance, 
could be calculated if K, and K, were evaluated. By burning various 
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substances in his lamp and measuring flame-heights and radii at their 
smoke points by means of a travelling microscope, he evaluated K, and K, 
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' and calculated A for several homologous series. From these calculated 
values of A the “‘ tendency to smoke ”’ was obtained from the equation : 


S, = 320/h 
His results are shown in Fig. 4 


Tasie II. 
Hydrocarbon Mizxtures. 





Smoke point 
Substance. flame-height, 


cm. 








26. Benzole (sp. gr. 0-88) 
27. Pool motor spirit (sp. gr. 0 738, aromatics 25‘ 
28. Aromatic-free Pool ab ot spirit ; 
29. Pool motor spirit distilled . . 
. Pool motor spirit distilled and freed from aromatics . 
. Aromatic-free Pool motor spirit distilled over the same range as 29 
2. Petroleum ether 40--60° C. ; ‘ ‘ 
32a, Aromatic-free petroleum ether 60-80° C. 
4 heii ial B.P. Solvent No. 


SaISSHP —SSAeas 
er OIwoo8sa ak wre Oo 


(aromatic- free) 





One of the most interesting features of these curves is the fact that the 
tendency to smoke of the hydrocarbon series investigated decreased with 
increase in the number of carbon atoms or boiling point, with the exception 
of the paraffins, which show an increased tendency to smoke with increase 
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in molecular weight. In the other series examined the increase in molecular 
weight is brought about largely by the addition of paraffinic side-chains, 
and the influence of the nucleus is thus progressively offset by the increased 
paraffinicity. 


Hydrocarbon Mixtures. 


The smoke points of commercial hydrocarbon mixtures are listed jp 
Table II. 

Pool motor spirit (27) is appreciably less smoky than benzole (26), and 
the removal of its aromatics produces further marked improvement, 
Results 29, 30, and 31 indicate that, if any high-boiling-point polymers 
are formed during the acid treatment, their effect on the tendency to smoke 
of the spirit is negligible. The 40-60° C petroleum ether (32) produces 
little smoke, but the aromatic-free 60-80° C fraction (32a) is rather more 
smoky. 

Of the Special B.P. Solvents, No. 1 is the least smoky (33), its flame. 
height being close to that of aromatic-free Pool motor spirit, while result 
33a shows that the removal of its aromatics (although only 3 per cent) 
produces an appreciable increase in flame-height, The flame-heights of 
Nos. 2 and 4 are very similar, with No. 5 a little more smoky. All of these 
special solvents have a greater flame-height than Pool motor spirit. 


OxYGEN-CONTAINING COMPOUNDS. 


(a) Alcohols. 
TABLE IIT. 


Alcohols. 








Oxygen in 
} Sommpound | sescke 
Boiling centage of point 
Substance. the total —- 
required for —~ 

complete = 
combustion. 








. Methyl alcohol Blue flame 
. Ethyl] alcohol 

. n-Propyl alcohol . 

. Isopropyl! alcohol 

. n-Butyl] alcohol 

. Isobutyl alcohol 

. Secondary butyl alcohol 

. Tertiary butyl alcohol . 

. n-Amy] alcohol 

. Isobutyl carbinol . : : - | 
. Secondary butyl carbinol ; 
. Methyl propyl carbinol (secondary amyl) | 
. Dimethyl-ethyl-carbinol ow amyl) . 

. n-Hexyl alcohol . ° 
. 1: 3-Dimethy] butyl alcohol . 

. 2-Ethyl butyl alcohol 

. Allyl alcohol , 

. Benzyl alcohol 

. cycloHexanol 

. Cresylic acid 
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With the straight-chain primary alcohols the tendency to smoke increases 
with increasing chain length, the flame-heights of the lower members 
showing a straight-line relationship with their oxygen content expressed as 
a percentage of the total required for complete combustion (Fig. 5). 





40 | 
ETHYL ALCOHOL 


4 





é 


N-PROPYL ALCOHOL ‘s 





YS w-puTvL ALCOHOL 


x 
° 





FLAME HEIGHT - CM 


ZL. N-AMYL ALCOHOL 
N-HEXYL ALCOHOL 


rd 


5 10 1s 
OXYGEN IN COMPOUND AS PERCENTAGE 
OF TOTAL REQUIRED FOR COMPLETE COMBUSTION 


Fie. 5. 














Further, for these compounds the inverse of flame-height shows a straight- 
line relationship with the boiling point (Fig. 6). The secondary alcohols 
are in all cases more smoky than the corresponding primary compounds, 
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which indicates that the shifting of the -OH group from the end of the 

chain to some other position increases the tendency to smoke. As in the 

case of the hydrocarbons, the branching of the chain increases the smoke, 

and it is interesting to note that a primary alcohol with a single-carbon 
ZZ ° 
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branch has roughly the same flame-height as the corresponding secondary 
alcohol ; results for isobutyl carbinol (46) and secondary butyl carbinol (47) 
suggest that the position of the branch relative to the primary -OH group 
has little or no effect on the smoking properties. 
(b) Glycol and Glycerol Compounds. 
TABLE IV. 
Glycol and Glycerol Compounds. 








Oxygen in 
compound | §moke 
Boiling as @ per- int 
Substance. point, centage of me 
*¢. the total | height, 
required for cm. 
complete 
combustion. 
57. Ethylene glycol . : ; ; P 197 28-6 Semi- 
luminous 
58. Cellosolve (ethylene glycol mono-ethyl 
ether) ‘ ‘ ; ‘ ‘ : 127 15-4 38-5 
59. Cellosolve acetate : ; 3 ; 156 16-7 28-2 
60. Ethylene glycol diacetate _.. : ' 188 23-5 29-2 
61. Triacetin . ; ; . : : 258 24-0 26-1 

















The importance of the C/H ratio in these compounds is stressed by a 
comparison of n-propyl alcohol (39) and allyl alcohol (53), since allyl 
alcohol, although having a higher oxygen content, is far more smoky than 
n-propyl alcohol. Also cyclohexanol has a greater tendency to smoke than 


TABLE V. 
Nitro-compounds. 

















Oxygen in 
Boilin Compound ¢ th Smoke point 
oili rcentage of the ; 
Substance. point, oo. ere required flame-height, 
for complete =, 
combustion. I 
62. Nitro-methane . . : 101 36-4 | Semi-luminous | 
63. Nitro-ethane : : , 114 : 23-5 44 (approx.) 
64. 1-Nitro-propane . .. 132 N_>NO 17-4 | 36-6 
65. 2-Nitro-propane . : ; 120 2) 17-4 20-3 ( 
66. Tertiary-nitro-isobutane : 128 13-8 7-9 
67. Nitrated Pool motor spirit . _— _ 5-2 
68. Nitrated aviation spirit . - = . q —~ 4:3 
69. Nitrated spirit 100 flat . -llpinte' agree a 4-5 | 
70. Nitrated spirit 90 flat { 8g | 





n-hexyl alcohol, although in this case ring formation may be an additional 
effect. Benzyl alcohol and cresylic acid, being essentially aromatic, are 
very smoky, in spite of their oxygen content. 

When the flame-height of cellosolve (58) is plotted against oxygen content 








dary 
| (47) 
Troup 
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in Fig. 8, the point lies close to the normal primary alcohol line, suggesting 
that the introduction of an oxygen atom into the carbon chain does not 
greatly modify the flame-height/oxygen content relationship. A com- 
parison of cellosolve (58) with cellosolve acetate (59) indicates that the 
presence of the acetate group increases the tendency to smoke and this 
receives some support from the comparison between ethylene glycol 
diacetate (60) and triacetin (61) (cf. Fig. 8). 

(c) Nitro Compounds.—The flame-heights of several pure nitro-com- 
pounds and of the products obtained by the nitration of some commercial 
hydrocarbon mixtures are tabulated in Table V. 

As in the case of aliphatic alcohols, the tertiary nitro-compound is 
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relatively smoky in spite of its high oxygen content. Also 2-nitro-propane 
(secondary) is more smoky than 1-nitro-propane (primary). 

Acetone and methyl ethyl ketone appear to have very similar smoking 
properties, although they differ in oxygen content. The flame-height of 
diacetone alcohol is much higher than might be expected, since it contains a 
tertiary alcohol grouping. Also the flame-height of cyclohexanone (which . 
is very close to that of cyclohexanol—9-8 cm) suggests that the substitution 
of O for H, in cyclohexane (7-7 cm) reduces the smoke. Phenyl-n-butyl 
ketone, although it contains a n-butyl chain, shows that the benzene ring 
is still strongly influencing the smoking tendency. 

The graph of smoke point flame-height against oxygen content for the 
straight-chain esters (Fig. 7) reveals an interesting relationship. For each 
set of isomers (members having equal oxygen content and similar boiling 
points in this case) there is a decrease in flame-height in the order butyrate, 
propionate, acetate, formate. Thus for each oxygen content the flame- 
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height at which smoking begins, increases with the chain length attached 
directly to the carboxylic carbon atom. Results for n- and iso-buty| 
acetate (86) and (87), for n- and iso-butyl propionate (94) and (95), and for 
m- and iso-butyl butyrate (100) and (101) indicate that, as with other 
compounds, the branching of a carbon chain tends to increase the smoke 
from these esters. Again results for ethyl formate and oxylate (80) and 
(103), and for ethyl acetate and malonate (84) and (104) suggest that the 
ethyl esters of dibasic acids have a greater tendency to smoke than the 
ethyl esters of mono-basic acids of similar oxygen content. A comparison 
of ethyl carbonate (110) and ethyl lactate (105) gives further evidence that 
the introduction of an oxygen atom into a carbon chain does not affect the 
flame-height /oxygen content relationship. 

Ethyl benzoate (108), although having a moderate oxygen content, is 
still very smoky, as would be expected. 

(d) Ketones. See Table VI. 

(e) Esters. See Table VIII. 

(f) Carboxylic Acids.—Tests were made on glacial acetic acid, oleic acid, 
and naphthenic acid, but, owing to the unsteady nature of the flame, 
reliable measurements were not possible: the first named appeared to 
have a high smoke point and the latter two, very low smoke points. 


Discussion. 


It is apparent from the above results that the tendency to smoke of an 
orgatfic compound on burning freely in air is intimately related to its 
molecular structure. Certain general relationships between smoking 
tendency and molecular structure may be expressed as follows. 


Hydrocarbons. 


Of all the hydrocarbons examined, the normal paraffins are the least 
smoky. With these compounds the tendency to smoke increases with 
increasing chain length. Although the normal paraffins have the lowest 
C/H ratios of the hydrocarbons, it is apparent that other factors also have 
an important bearing on the smoking tendency. Thus branched-chain 
paraffins are more smoky than their normal isomers, although they have the 
same C/H ratios. Indeed, the more highly branched the carbon chain, the 
more smoky does the paraffin become, while the relative positions of the 
branches appear to have little effect. 

Although the olefines examined were probably mixtures, the results 
suggest that unsaturation in a hydrocarbon produces an appreciable 
increase in the tendency to smoke. It seems probable, however, that an 
increase in chain length of an olefine results in a decreased smoking tendency, 
since the effect of the double bond is progressively offset by increased 
paraffinicity. 

The tendency to smoke of the naphthenes, although relatively high in 
comparison with the normal paraffins, is approximately equal to that of 
those paraffins having two or three branches in the carbon chain. There is 
insufficient evidence to show the effect of the number of carbon atoms in 
the naphthene ring on the smoking tendency, but a comparison of cyclo- 
pentane and cyclohexane suggests that the effect is not pronounced. As 
with other hydrocarbons, the introduction of one double bond in a napthene 
produces a marked increase in the tendency to smoke. 
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TasLe VI. 
Ketones. 
Oxygen, in 
compound as @ , 
An Smoke point 
Substance. point, #0 pore pede flame-height, 
for complete =. 
combustion. 
7. Acetone 56 11-1 26-3 
72. Methyl-ethyl ketone 80 8-3 26-2 
73. Di-acetone alcohol 164 11-1 26-4 
14. cycloHexanone 155 5-9 10-0 
15. } n-butyl ketone 246 3-4 1-2 
Tasie VIII. 
Esters. 
Oxygen in 
Compound 
as & per- Smoke 
Substance. Boiling centage of 4 
point, ° C the total — 
required for | height, 
complete = 
combustion 
80. Ethyl formate / 54 22-2 38-5 
81. n-Propyl formate . 81 16-7 15-6 
82. n-Butyl formate . 107 13-3 11-6 
83. Methyl acetate 57 22-2 About 45 
84. Ethyl acetate 77 16-7 28-6 
85. n-Propyl acetate 101 13-3 13-0 
86. n-Butyl acetate 126 11-1 10-5 
87. Isobutyl acetate 116 11-1 9-2 
88. Amyl acetate : 131-135 9-5 10-0 
89. cycloHexyl acetate 175 8-7 9-2 
90. Octyl acetate (2-ethyl- hexyl acetate) 196 6-7 9-6 
91. — acetoacetate 181 17-6 21-6 
92 re propionate 79 16-7 36-0 
93. Ethy propionate 4 98 13-3 24-5 
94. n-Butyl propionate 146 6 9-5 13-7 
95. Isobutyl propionate 137 9-5 10-3 
96. Amyl propionate 156 8-3 10-5 
97. Meth I eatgreate 102 13-3 27-5 
98. Ethyl butyrate 120 11-1 22-9 
99. Propyl butyrate 143 9-5 14-4 
100. n-Butyl butyrate 165 8-3 13-3 
101. Isobutyl butyrate . 157 8-3 11-3 
102. Ethyl laurate » 269 4-8 18-4 
103. Ethyl oxalate 186 23-5 29-0 
104. Ethyl malonate 198 20-0 25-7 
105. Ethyl lactate 154 20-0 32-1 
106. Butyl lactate 187 14-3 19-0 
107. Amy] lactate 198 12-5 15-4 
108. Ethyl benzoate 213 8-7 2-5 
109. Ethyl nitrate 87 35-3 Semi- 
luminous 
110. Ethyl carbonate 126 20-0 33-0 
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Benzene, having three double bonds in its ring and a correspondingly high 
C/H ratio, is very smoky, and the building up of this effect may be traced 
from n-hexane through cyclohexane to cyclohexene and benzene. Further, 
the results obtained with this lamp indicate that aliphatic side-chains on 
the benzene ring have little effect on its tendency to smoke. 

In the naphthalene series the effect of unsaturation (or C/H ratio) is 
shown, since decalin (C/H ratio 0-56) is appreciably less smoky than 
tetralin (C/H ratio 0-83). 


The Compactness of a Hydrocarbon Molecule and the Tendency to Smoke. 


It seems likely that when a hydrocarbon is decomposed in a flame, the 
more compact the molecule, the greater will be the initial size of the carbon 
particles. This in its turn will probably produce a-greater tendency to 
smoke. Thus a branched-chain paraffin, having a more compact molecule, 
is more smoky than its straight chain isomer, although they have the same 
C/H ratio. Ring formation also produces a compact molecule, and these 
substances, even when saturated, are more smoky than the corresponding 
straight chain compounds. 








Oxygen-Containing Compounds. 
The effect of the structure of oxygen-containing organic compounds on 
the tendency to smoke on burning may be resolved into three factors : 


(1) The structure of the hydrocarbon portion of the molecule ; 

(2) The structure of the oxygen-containing group or groups ; 

(3) The oxygen content as a percentage of that required for complete 
combustion. 


The structure of the hydrocarbon portion of the molecule is of primary 
importance, and largely determines the tendency to smoke of the compound. 
It would appear that the relationships between the structure of the hydro- 
carbon portion and the smoking tendency of the substance are largely the 
same as for the hydrocarbon compounds. In general, increased oxygen 
content results in a decreased tendency to smoke, but this effect may be 

‘ offset by the influence of the hydrocarbon portion. ° 

With the straight-chain primary alcohols the tendency to smoke increases 
with increasing chain length, the flame-heights of the lower members 
showing a straight-line relationship with their oxygen contents expressed 
as a percentage of the total required for complete combustion. Further, for 
these compounds the inverse of flame-height shows a straight-line relation- 
ship with the boiling points. Tertiary alcohols are more smoky than the t 
corresponding secondary alcohols, which are in turn more smoky than 
the primary compounds. 

The importance of the hydrocarbon structure in alcohols is apparent in 
several instances. Firstly an alcohol with a branched chain is more smoky 
than the corresponding straight-chain compound, although they have the 
same oxygen content. Also an alcohol having an unsaturated carbon 
chain is, in spite of its higher oxygen content, far more smoky than the 
corresponding saturated alcohol. Further, a saturated cyclic alcohol has 
a greater tendency to smoke than the straight-chain compound, and an 
aromatic alcohol shows the pronounced influence of the aromatic ring in 

spite of its oxygen content. 
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With the straight-chain aliphatic esters the tendency to smoke is 
influenced by the carbon chain length attached to the carboxylic carbon 
atom and to the hydroxyl oxygen atom respectively. For each set of 
isomers (members having equal oxygen content and similar boiling points 
in this case) there is a variation in the smoke point flame-height, and the 
order of decreasing tendency to smoke is: formate, acetate, propionate, 
butyrate. Thus for equal oxygen content the flame-height at which 
smoking begins increases with the chain length attached directly to the 
carboxyl carbon atom. 

As with the other types of compounds, the branching of a carbon chain 
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tends to increase the smoke from these esters. Also the results suggest 
that the ethyl esters of dibasic acids are more smoky than the ethyl esters 
of monobasic acids of similar oxygen content. Aromatic esters are again 
very smoky. 

It would appear from the few aliphatic nitro-compounds examined that 
the tertiary compound is more smoky than the corresponding secondary 
compound, which is in turn more smoky than the primary compound. 

It is interesting to note that with several substances the introduction of 
one or two oxygen atoms into the carbon chain appears to have little or no 
effect on the flame-height /oxygen content relationship. 

The relative smokiness of the various types of organic compounds may 


be observed from Fig. 8. The introduction of an oxygen atom into n- 
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pentane (1) to form n-amyl alcohol (45) produces a marked increase i 
flame-height. If, however, the flame-heights of n-hexane and n-hepty] 
alcohol are obtained by extrapolation, the results suggest that the increase 
from n-hexane to n-hexyl alcohol (50) is less marked than for the aboyg 
5-carbon compounds and that the conversion of n-heptane (4) into n-hepty 
alcohol produces no increase in flame-height. It is interesting to note that, 
due to their structures, a number of compounds containing fairly high | 
percentages of oxygen are more smoky than hydrocarbons, such as n 
pentane (1) and n-heptane (4). Examples of such substances are : 


Unsaturated aliphatic alcohol-allyl alcohol (33) ; 
Nitro-paraffin-tertiary nitro-isobutane (66) ; 

Aliphatic esters—n-butyl formate (82) and n-propyl acetate (85) ; 
Aromatic esters—ethyl benzoate (108). E 

However, n-pentane and n-heptane (the best of the hydrocarbons 
examined) are far more smoky than the best of the oxygen- containing ~ 
compounds examined. Of the latter, the normal primary alcohols have the | 5 
greatest flame heights on an equal oxygen content basis. 

Although methyl ethyl ketone (72) lies above the alcohol line, acetoniil 
(71), having the highest oxygen content (11-1 per cent) of this ketone series, ” 
lies below the line. Nitro-ethane (63) and 1-nitro-propane (64) are more — 
smoky than the normal primary alcohols for equal oxygen content. The 
oxygen content of a mono-nitro paraffin is, however, much higher than that 
of the alcohol having the same number of carbon atoms. Hence 1-nitro- 
propane (64) hag a much greater flame-height than n-propyl alcohol (39), 
and, similarly (only to a lesser degree due to the different slope of the 
lines), nitro-ethane (63) is less smoky than ethyl alcohol (38), while. by 
extrapolation (assuming for this purpose a straight-line relationship for the 
nitro-paraffins) the flame-height of nitro-methane (oxygen 36-4 per cent of 
total required) would be roughly equal to that of methyl alcohol (oxygen 
25-0 per cent of total required). 

For any particular set of isomeric aliphatic esters, the methyl ester of the 
highest acid has the greatest flame-height. However, the results suggest 
that of these methyl esters, methyl formate is the least smoky, due to its 
high oxygen content. 

At equal oxygen content the general order for increasing tendency to 
smoke is : 


(1) Normal primary alcohols ; 

(2) Normal primary nitro-paraffins ; 
(3) Propionates ; 

(4) Acetates ; 

(5) Lactates ; 

(6) Formates ; 


although this order will vary slightly for different oxygen contents. 
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